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Abstract —We consider the problem of coordinated multi¬ 
cell downlink beamforming in massive multiple input multiple 
output (MIMO) systems consisting of N cells, N t antennas 
per base station (BS) and K user terminals (UTs) per cell. 
Specifically, we formulate a multi-cell beamforming algorithm 
for massive MIMO systems which requires limited amount of 
information exchange between the BSs. The design objective is 
to minimize the aggregate transmit power across all the BSs 
subject to satisfying the user signal to interference noise ratio 
(SINR) constraints. The algorithm requires the BSs to exchange 
parameters which can be computed solely based on the channel 
statistics rather than the instantaneous CSI. We make use of 
tools from random matrix theory to formulate the decentralized 
algorithm. We also characterize a lower bound on the set of target 
SINR values for which the decentralized multi-cell beamforming 
algorithm is feasible. We further show that the performance 
of our algorithm asymptotically matches the performance of 
the centralized algorithm with full CSI sharing. While the 
original result focuses on minimizing the aggregate transmit 
power across all the BSs, we formulate a heuristic extension of 
this algorithm to incorporate a practical constraint in multi-cell 
systems, namely the individual BS transmit power constraints. 
Finally, we investigate the impact of imperfect CSI and pilot 
contamination effect on the performance of the decentralized 
algorithm, and propose a heuristic extension of the algorithm 
to accommodate these issues. Simulation results illustrate that 
our algorithm closely satisfies the target SINR constraints and 
achieves minimum power in the regime of massive MIMO 
systems. In addition, it also provides substantial power savings 
as compared to zero-forcing beamforming when the number of 
antennas per BS is of the same orders of magnitude as the 
number of UTs per cell. 

Index Terms —Massive MIMO, coordinated beamforming, de¬ 
centralized design, random matrix theory. 

I. Introduction 

Massive multiple input multiple output (MIMO) has been 
identified as an essential ingredient in the design of next gener¬ 
ation cellular systems, as it provides substantial improvement 
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in both spectral and energy efficiency m. It refers to the idea 
of scaling up the number of antennas on the base station (BS) 
to a few hundreds, serving many tens of user terminals (UTs) 
on the same resource block. The basic idea is to exploit large 
number of antennas to achieve greater spatial resolution and 
array gain, resulting in a higher throughput and greater energy 
efficiency. 

The idea of massive MIMO was first proposed in the 
seminal work of E). The main finding of the paper was 
that as the number of antennas on the BS grows without 
bound, the effects of fast fading and interference vanish, 
and the system performance is ultimately limited by only 
pilot contamination 0. The other attractive feature is that 
simple signal processing techniques at the BS, such as the 
use of eigen beamforming and matched filter were optimal 
under this setting. Subsequently, massive MIMO systems was 
also studied from an energy efficiency point of view and 
shown to achieve dramatic improvement in this regard 0 . The 
impact of channel estimation, pilot contamination, and antenna 
correlation in massive MIMO systems was investigated in 0. 
In particular, it was concluded that sophisticated beamforming 
techniques such as regularized zero forcing (RZF) outperform 
eigen beamforming in a massive MIMO setting when the 
number of antennas is of the same orders of magnitude as 
the number of users. Reference 0 addressed the question 
of how to select the system parameters in a massive MIMO 
system (number of antennas per BS, number of users, transmit 
power etc.) to maximize the energy efficiency. The work in 
0 proposes a hierarchical interference mitigation scheme in 
massive MIMO systems based on a two level precoding with 
the objective of maximizing the system utility. Although all 
these studies point to impressive gains in massive MIMO both 
in terms of spectral efficiency and energy efficiency, construct¬ 
ing such large dimensional arrays can result in significant 
additional hardware cost of the analog front ends. Moreover, 
extra physical dimensions are required in order reduce the 
mutual coupling between the antenna elements. 

Subsequently, it was proposed that in a multi-cellular envi¬ 
ronment, the gains obtained by massive MIMO systems can be 
replicated by using much lesser number of antennas if BSs are 
allowed to cooperate with each other. In this context, 0 pro¬ 
posed a TDD architecture based network MIMO like system 
with BS cooperation and zero-forcing (ZF) beamforming, and 
showed that massive MIMO performance can be achieved with 
one order of magnitude fewer antennas per active user per cell. 
These results motivate us to consider multi-cell cooperation 
in massive MIMO systems. However, since massive MIMO 
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systems are inherently large, enabling BS cooperation in a 
multi-cellular environment requires tremendous amount of 
information exchange between them. 

In order to address this issue, in this work, we propose 
an optimal decentralized multi-cell beamforming algorithm 
for massive MIMO systems that requires limited amount of 
information exchange between the BSs. We primarily focus 
on the so called coordinated beamforming, in which BSs 
formulate their beamforming vectors taking into account the 
interference they cause to the neighboring cells This is 
accomplished by exchanging the CSI information between 
them. However, unlike network MIMO systems, no user data 
exchange takes place between the BSs. The design objective 
considered is to minimize the aggregate transmit power across 
all the BSs subject to satisfying the user signal to noise 
ratio (SINR) constraints. Reference 0 provides an optimal 
centralized algorithm to solve this problem. However, the 
centralized solution demands high computational ability, and 
exchange of the fast fading CSI co-efficient between the 
BSs. Such an algorithm requires high capacity backhaul links, 
especially when implemented in a massive MIMO setting. 

In order to overcome the heavy backhaul requirement, we 
propose in this work a decentralized approach to compute the 
multi-cell beamforming vectors. In our algorithm, the BSs 
must exchange parameters at the time scale of slow fading 
coefficients rather than the instantaneous channel realizations 
(fast fading coefficients). We use tools from random matrix 
theory (RMT) to formulate our algorithm. 

The multi-cell beamforming strategy involving exchange of 
parameters based on channel statistics was first proposed in 
m using tools from RMT. It was shown with the help of sim¬ 
ulations that such an algorithm performs well for large system 
dimensions. Similar ideas for multi-cell beamforming were 
subsequently proposed in CD, 02 and theoretical arguments 
for the asymptotic optimality of this algorithm were provided 
in the special case of a two cell Wyner model, with symmetric 
SINR constraints for all the UTs. However, the analyses in 
CD, 02 rely on obtaining closed form expressions for the 
system parameters (such as the uplink and downlink power), 
and are not extendable to more practical channel models. In 
this work, we provide a comprehensive design of the RMT 
based decentralized beamforming under a massive MIMO 
multi-cell setting with a distance based pathloss model. Fur¬ 
ther, we provide arguments for asymptotic optimality of this 
algorithm. Specifically, our main contributions are as follows: 

For the original problem of minimizing the aggregate trans¬ 
mit power across all the BSs subject to UT SINR constraints, 
we present the following results: 

• We propose a reduced overhead beamforming algorithm 
(ROBF) in a massive MIMO multi-cell setting. In this 
algorithm, the BSs require the knowledge of local CSI 
(of the UTs they are serving and also the UTs present in 
the other cells). In addition, the BSs have to compute 
parameters that depend only on the channel statistics, 
which they exchange between them to compute the 
beamforming vectors. 

• Using a large system analysis, we provide closed form 
expression for the lower bound on the set of target SINR, 


for which the decentralized multi-cell beamforming algo¬ 
rithm is feasible. 

• We prove that when the dimensions of the system become 
large, the achieved SINR in the uplink and downlink 
by the ROBF algorithm exactly match the target SINR. 
Moreover, we also prove that when the dimensions of the 
system become large, the performance of our algorithm in 
terms of uplink and downlink transmit powers perfectly 
match that of the optimal algorithm proposed in 0. 

With this algorithm as reference, we present heuristic ex¬ 
tensions to incorporate two practical constraints in MIMO 
multi-cell systems (1) individual BS transmit power constraints 
(2) impact of imperfect CSI and pilot contamination. The 
contributions are as follows: 

• We formulate a heuristic extension of the decentral¬ 
ized multi-cell beamforming algorithm to incorporate the 
individual BS transmit power constraints, and present 
numerical results to show the convergence as well as the 
performance of this algorithm. 

• Finally, we investigate the impact of CSI estimation and 
pilot contamination on the performance of the ROBF 
algorithm and propose a heuristic adaptation of the ROBF 
algorithm that can provide better performance in the 
presence of pilot contamination. 

In addition, our work contains several novel ideas of combin¬ 
ing RMT results with optimization theory and uplink-downlink 
duality in MIMO systems, which are of independent interest. 

It is worth noting that RMT results have been used exten¬ 
sively to assess the performance of linear beamforming strate¬ 
gies in multi-user MIMO systems lH3l . [fUfl l. [3), J8]|. However, 
all these works consider “pre-defined” transmit strategies such 
as eigen beamforming, zero forcing (ZF), regularized zero 
forcing (RZF) etc., and use RMT as a tool for performance 
analysis of these schemes in the large dimensional regime. In 
contrast, in this work, we use RMT results for the purpose 
of system design (as opposed to performance analysis), i.e., 
design of optimal beamforming vectors in MIMO multi-cell 
system^] Such an approach is novel and is facilitated by 
combining RMT results with optimization theory. Moreover, 
optimizing the system performance imposes additional techni¬ 
cal difficulties in applying RMT results. For e.g., implementing 
a power control algorithm (both in the uplink and downlink) 
implies that the transmit powers explicitly depend on the 
channels (and hence the randomness associated with the chan¬ 
nel realization). Such dependency of transmit powers on the 
channel realization makes it unsuitable to apply RMT results. 
Our approach in this work is to first propose an algorithm 
that depends only on the second order statistics of the channel 
vectors (the path-loss in our case). Then, we apply such an 
algorithm to the original system set-up, and prove that this 
algorithm is optimal in the large system domain. 

Although the theoretical results prove the optimality of 
ROBF algorithm in the asymptotic regime, we provide nu¬ 
merical results to show that the performance of ROBF algo¬ 
rithm closely matches that of the centralized algorithm for 

'in an unrelated context, RMT results have also been used in the context 
of system design in on and COD 
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moderate system dimensions (when the number of antennas 
are comparable to the number of UTs per cell), both in 
terms of satisfying the user SINR targets and minimizing the 
downlink transmission power. Moreover, these results indicate 
that ROBF algorithm provides substantial trasnmit power 
savings as compared to other beamforming strategies such as 
zero-forcing n the regime where the number of antennas is 
comparable to the number of UTs. 

Finally, we remark that all the analysis in this work is 
performed assuming independent and identical (i.i.d.) channel 
vectors and time division duplex (TDD) mode of operation. 
The performance of massive MIMO systems with correlated 
channel models have been studied in prior works such as 
0 . Recently, there has also been an interest in exploring 
massive MIMO systems with frequency division duplex (FDD) 
mode of operation (7), ifTTI . fTSI . As this work is a first step 
towards exploring power control and the design of optimal 
beamforming in massive MIMO systems, in order to keep 
the analysis simple, we restrict our attention to i.i.d. channels 
and TDD mode of operation. The extension to the case of 
correlation channels and FDD mode of operation will be a 
topic of future research. 


II. System Model and Algorithm Description 
A. System Model 

We consider the problem of multi-cell beamforming across 
N cells and I\ UTs per cell where each BS is equipped with 
Nt antennas and each UT has a single antenna. Each BS 
serves only the UTs in its cell. Let h€ C■ :V ' denote the 
channel from the BS i to the fc-th UT in cell j. We consider 
reciprocity between the uplink and downlink channels, and 
hence the TDD mode of operation, as it is the preferred 
mode of operation in massive MIMO systems 0. We assume 
that the elements of the channel vector are independent and 
identically distributed (i.i.d.) with Gaussian distribution, i.e., 
h ij t k ~ CM(0, <Jij tk I Nt ), the variance (Jij.k of the channel 
depends upon the path loss model between BS i and UT(j, k). 
Recent works on channel measurements indicate that i.i.d. 
assumption is a reasonable model for massive MIMO arrays 
m. We assume that the BSs have perfect CSI of the downlink 
channels to all the users in the system (i.e., h iiTlj fc, Vn, k). Let 
w i j £ C Nt denote the transmit downlink beamforming vector 
for the 7 -th UT in cell i. Likewise, let A™, denote the received 
SINR for the jth UT in cell i and jij the corresponding target 
SINR. The received signal £ C for the jth UT in cell i, 
is given by 


The rest of the paper is organized as follows. We provide 
the system model and describe our reduced overhead multicell 
beamforming in section [D] In section [III} we provide the 
asymptotic analysis of the reduced overhead algorithm formu¬ 
lated in the previous section. In Section [VJ we investigate the 
impact of imperfect CSI and pilot contamination on the perfor¬ 
mance of the ROBE algorithm. We summarize the simulation 
results in section [VI] Linally, we provide concluding remarks 
in section VII Appendix A provides some relevant results 


from RMT which will be used in formulating our algorithm. 
Appendices B, C, D, E, L, G and H provide the proofs of some 
of the results stated in the paper. 


Notations: Throughout this work, we use boldface lower¬ 
case and uppercase letters to designate column vectors and 
matrices, respectively. Lor a matrix X, X(p, q) denotes the 
( p,q ) entry of X. X T , X ff , tr(X), ||X||, and p(X) denote 
the transpose, the complex conjugate transpose, the trace, 
the spectral norm and the spectral radius of the matrix X, 
respectively. Lor two matrices X and Y, the notation X < Y 
denotes element wise inequality (X(p, q) < Y (p,q) Vp,q) 
and similarly for vectors. We denote an identity matrix of 
size M as I m and diag(aq,..., Xm) is a diagonal matrix of 
size M with the elements 27 on its main diagonal. We use 
x ~ CM {m, R) to state that the vector x has a complex 
Gaussian distribution with mean m and covariance matrix R. 
We use the notation —> to denote almost sure convergence. 
Let ay and 6y denote a pair of infinite sequences. We write 
a,N x &yr, iff ayr — feyr 0. We denote the expectation of 
a random variable by the notation E [.] Let N t ,K £ N + , we 
use the notation N t , K —► 00 to denote the following condition 
on N t and K , 0 < liminf 00 77 < limsup^.,.^ ^ < 00 . 
Linally, the notation ( x) + is used to denote max(a;, 0). 


Vi,: 


= h^-w i 


T 


E 

(n,k)^(i,j) 


h 


,k x n,k T 


where Xij £ C represents the information signal for the j- 
th user in cell i and z i: j ~ CM (0. N 0 ) is the corresponding 
additive white Gaussian complex noise. Under this model, the 
achieved SINR in downlink for the UT ij is given bvj^] 


A = — 

l ’ J E 






H 


.12 


- No 


(1) 


The denominator terms of {!]) represent the intra-cell inter¬ 
ference, inter-cell interference and the noise (in order as they 
appear). The downlink sum power minimization problem can 
be formulated as the following optimization problem given by 


min 

W ij 

S.t. 


E 

ij 


A DL > 'y- ■ 
1Y z,j — 7 




( 2 ) 


B. Algorithm Design 

As show in l20l . the optimization problem 0 can be 
reformulated as a second order conic programming (SOCP) 
problem and, strong duality holds for this problem. Lollowing 
the approach of a, we solve this problem using duality 
theory. Accordingly, we introduce the Lagrange multiplier 
-4B- associated with the downlink SINR constraints. The 

JAlt 

2 Note that the downlink SINR expression in {T| assumes that the channels 
are known perfectly at the UTs. However, the main focus of this work is in 
the massive MIMO regime, i.e., Nt and K being very large. Fortunately, as 
shown in 0, under this regime, the UTs only need to have the knowledge 
of the average effective channels. 
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Lagrangian is given by 


L( w , A) = £ w£-w itj - ^ ir 


A ij r | 

iV* L 7,^- 


1,3 

^ [ |w n ^h n (V[ 

(n,k)jt(i,j) 


(3) 


Note that the Lagrange multiplier is scaled by the factor N t 
in order to ensure that the sum power in the system is finite, 
when the dimensions of the system grow large (in terms of the 
number of antennas on the BS and number of UTs). In fact, 
the Lagrange multipliers can be interpreted as the dual 
uplink powers in the formulation of the dual uplink problem 
obtained in the following manner. 

Rearranging ([3|, we obtain 


L( w, A) 


E AjjiVo 

N t 


hj 


+ E w / // . 

h3 


(4) 


where the matrix is given by 



= I 


'*,3 






E 

(n,k)^,(i,j) 


A n,k u ,H 

Nt ai , n ’ kn i,n,k- 


(5) 


The dual uplink problem corresponding to the optimization in 
0 is formulated as 


min 

Ai.j, Vi,j 

s.t. 



Vi, j 


( 6 ) 


where the left hand side of the constraint equation represents 
the uplink SINR given by 


A ul = 

l,3 


Ai,j 

N t 




E A n,k 

(n,fc)#(«,j) N t 




i,n,k | 
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where w l;J denotes the corresponding uplink receive filter. 

We now provide a brief description of the beamforming 
algorithm presented in f9|. Before introducing the algorithm, 
we define the following matrices. 


N t xK 


Hz,n — [hi,7i,l, ■ • • , hz,n,lv] tz 

Hj = [H ijl5 ..., H^jv] G C NtXNK 


Aj — [^2,1) * • • ■) ^2,.ft”] ^ 

A = diag [Ai,...,Ajv] e 


iKxl 


xNKxNK 


We also define the matrix S A = G C N±xNt 

L l\t L 


Algorithm 1 (Centralized Algorithm - CBF). Perform the 
following steps. 


Starting from any initial A° • > 0 Vi, j the uplink power 


allocation is given by Xij = lim^oo A* j, where 

1 


AJ+ 1 = 


v,n • • Tv.) 'ic,. y 


Vi, j 


(7) 


where = ^HiA*Hf and A ( 

diag [A*,..., A'y] . 

The optimal receive uplink receive filter is given by 

1 <E^WhL*+V„l)" 1 h,„ J . 


VN t ' 


( 8 ) 


The optimal transmit downlink beamforming vectors are 
given by w,j = J ^Wij, where Sij is given as 


d = F _1 liVn. 


Here, 


Kxl 


&i ~ • • • , € 

8 = [S 1 ,...,S n ]g C NKx1 
le [l,...,i] T e R NKx1 

and the elements of the matrix F € C NK x NK and the 
submatrix F l,J G C KxK are given by, 


F 1,1 ... F 1,JV 


F = 


(9) 


pJV,l 


... F 


N,N 


F i,n A f ^Iw^.hi, ijl 2 , n = i, k = j 
^\™n,k h n,i,j\ 2 , 


■ j,k 


We remark that the scaling of \fN~t in the expression for 
the uplink receive filter <|5}, and in the definition of the scaling 
factor 5ij ensure that the total power in the system is finite, 
when the dimensions of the system grow large. 

As mentioned before, the solution provided in m cannot 
be implemented in a distributed manner. The computation of 
dual uplink power (A, 7 ) and the scaling factors (<5,.y) requires 
a central station which has the global CS1 knowledge. In what 
follows, we overcome this problem. 

We now formulate our reduced overhead beamforming al¬ 
gorithm. The main idea behind this algorithm is that under the 
massive MIMO regime (i.e. when N t and I\ become large), 
the parameters in |7]i and ([TO]) can be approximated by their 
asymptotic equivalents using results from RMT. Moreover, 
the computation of these parameters will only depend on the 
second order statistics (path-loss), and not on the fast fading 
component of the channel vectors. However, note that RMT 
results are not directly applicable to this scenario. This is 
due to the fact that the computation of Xij in ([7]) explicitly 
depend on the channel vectors (RMT results require that the 
matrix A in 0 are independent of the channel matrices). This 
imposes additional technical difficulties in the application of 
RMT results in our scenario. Our approach in this work is to 
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first propose an algorithm that depends only on the second 
order statistics of the channel vectors (the path-loss in our 
case). Mathematically/theoretically, it is not ensured that it 
achieves the optimal solution. However, we apply such an 
algorithm to the original system set-up, and prove that this 
algorithm is optimal in the large system domain. 

We hereby represent the dual uplink power, the uplink and 
downlink beamforming vectors of the decentralized algorithm 
by the notation p-ij , v, :j and j ; respectively, which are the 
counterparts of A, j, w, ; and w, ; of the CBF algorithm. 


m' can be evaluated from rrii as 


-1 1 \ - ^ (&i,n,kUn,k'mi) 

1 Nt l^n=l zLfc=l (1+cr 


and the term 


Gi,n,k 


*L,ra,fc 


2 ' 


The vector p r = 

P = [Pii P 21 • • • i PnV £ 


(1 T ltn,k®i,n,ktTli) 

N 0 _ No 


»NKxl 


(17) 


(18) 


and 


Algorithm 2 (Reduced Overhead Beamforming algorithm - 
ROBF). Perform the following steps. 

• Starting from any initial jf -> 0 Vi, j the uplink power 


allocation is given by pij = lim^oo p\ •, where 


( 11 ) 


and m\ is evaluated as m\ = lim^oo m*’ p (initializing 
with any rh > 0,Vi) 


N K 






^i,n,kP n k 


t 1 


* „=1 fc=l 1 + a ^n,kP n:k m. 


+ 1 


( 12 ) 


The optimal receive uplink receive filter is given by 

(13) 


n,k 


The optimal transmit downlink beamforming vectors are 
given by Vjj = \ I jf-Vij. The scaling factor 5i t j is given 


as 


where 


5 = (i-rA)- 1 P , 


(14) 


Si = [Si,i,5 it 2,...,S <K ] T e M Xxl 

5 = [g 1 ,5 2 ,...,5 Ar ] T GM Arifxl 

7i,i 7i,*r 


7i = 


1 T 


7 = [7i> • • • i 7 jv] T 
T = diag(~f) 

and the matrix A G C ' NK xNK ^ e j^ ne( ^ as 


A = 


A 1 ’ 1 ... A 1 ^ 


A^ 1 _ A JVW 


where each submatrix A 1 ’-' G C KxK is given by 


(15) 


^ = 


0, n = i, k = j 

7 y 7 G n .%.jG n _ n ^ktti n ■ (n,k) 


(16) 


We now provide some remarks on this algorithm. 

Lemma 1. 77;c iterative algorithm © converges to a fixed 
point. 

Proof: The proof is provided in Appendix B. ■ 

In section III, we characterize the solution provided by this 
fixed point, and show that it is asymptotically optimal in the 
sense that when the dimensions of the system grow large, the 
achieved uplink and downlink SINR by this algorithm exactly 
match the target SINR, and the allocated uplink and downlink 
powers match the optimal solution. 

In the rest of the paper, we address the decentralized beam¬ 
forming algorithm by the acronym ROBF. We now discuss 
the practical advantages of the ROBF algorithm over the CBF 
algorithm. 


C. ROBF Algorithm - Signaling Overhead and Complexity 
Reduction 


We now provide a brief discussion on the signaling overhead 
and complexity reduction associated with the ROBF algorithm. 

Signaling Overhead: 

Recall that in the ROBF algorithm only the statistical infor¬ 
mation (path loss) must be exchanged between the BSs where 
as in the CBF algorithm, the fast fading co-efficients need 
to be exchanged. Let us denote the channel coherence time 
by t co h units, and the long term time constant by tlt units 
(over which the statistical properties of the channel change). 
In has been demonstrated in some prior works related to 
channel measurements, that in an urban setting the statistical 
information of the channel can be viewed as constant roughly 
for over 100 channel coherence intervals ED. 

We now characterize the signaling overhead for the two 
algorithms. For the CBF algorithm, exchanging full CSI would 
imply NN t K complex channel coefficients of 2NN t K real 
numbers every r co h units of time. This would amount to 
exchanging NN t K complex channel coefficients or 2NN t K 
real co-efficients. Therefore, the rate of information exchange 
for CBF algorithm would be 


f?CBF 


2 NN t K 

Loh 


real coefficients/sec. 


For the ROBF algorithm, NK real numbers (path loss in¬ 
formation assuming i.i.d. channel model) must be exchanged 
every tlt units of time, and the resulting rate of information 
exchange is given by 


7?robf 


NK 

r LT 


real coefficients/sec. 
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In order to get real feel of these numbers, assume TV = 3, N t = 
100, I< = 40, r LT = 22.6s, r coh = 180ms EJ. Therefore, 

/T( bp = 1-3 x 10 5 real coefficients/sec; 

TTrobf = 5.3 real coefficients/sec, 


Therefore, the ratio of two quantities can be given by 


T?cbf 

T?robf 


= 2 N t 



10 4 . 


However, for correlated channel model, T?robf 
NNtK real coefficients/sec. Therefore 


T?cbf _ 2 ( Tlt 
TTroBF \7coh 


100 . 


Finally, note that the exact number of bits to be exchanged 
depends on how these channel co-efficients are quantized 
(which is beyond the scope of this paper). 

Implementation Complexity: 

An exact characterization of the complexity associated with 
the ROBF algorithm is beyond the scope of this work. Nev¬ 
ertheless, we provide a brief analysis of the same. 
Computing the uplink power allocation 
Recall that implementing the iterations for the computation 
of the uplink power in the CBF algorithm requires matrix 
inversion operations, where as the ROBF algorithm only 
requires performing scalar operations. This tremendously re¬ 
duces the computational complexity with respect to the ROBF 
algorithm (the complexity of inverting a matrix is provided 
next). Moreover, in the CBF algorithm the uplink power 
must be evaluated for every channel realization (fast fading 
CSI), where as the ROBF algorithm requires parameters to 
be computed only once (at the time scale of changing of 
slow fading CSI). Therefore, there is a huge reduction in the 
computational complexity. 

Formulation of the downlink beamforming vectors 

Note that the downlink beamforming vector in the ROBF 
algorithm (and also the CBF algorithm) is in the form of a 
regularized zero forcing (RZF) beamforming, which requires 
computation of a matrix inverse of dimension N t x NK. 
Therefore, its computational complexity scales as N t (NK) 2 . 
This can be computationally demanding especially in a mas¬ 
sive MIMO setting. Fortunately, alternate schemes are being 
developed to implement RZF based on truncated polyno¬ 
mial expansion incur much less computational complexity 
as compared to matrix inversion, which are shown to have 
performance very close to RZF beamforming l22l . 

The price to pay for the reduction in information exchange 
between the BSs is that in the ROBF algorithm, the target 
SINR values are not met perfectly for every channel realiza¬ 
tion. In fact, the achieved SINR in the downlink fluctuates 
around the target SINR. However, we show through simula¬ 
tions that even for practical values of the system dimensions, 
the fluctuations of the achieved SINR around the target SINR 
value are small. In order to make sure that the target SINR 
requirements are satisfied, one could solve the optimization 
problem with ROBF algorithm by considering a higher value 
target SINR (than the actual desired one) in order to compen¬ 
sate for the fluctuations. 


Next, we show that the performance of ROBF algorithm 
perfectly matches the CBF algorithm when the number of 
antennas per BS and the number of UTs become large, i.e., 
in the regime of massive MIMO systems. We also provide 
simulation results to examine the performance of the ROBF 
algorithm. 


III. Algorithm Analysis 

In this section, we provide extensive analysis of the ROBF 
algorithm. In the rest of the paper, we use the phrase “large 
system” to refer to the regime when the number of antennas 
per BS and the number of UTs per-cell become large, i.e., 
N t ,K —)• oo while their ratio -jf tends to a finite constant 
0 < (3 < oo, as considered in some past works in this field 
0 . However, we mention that our results provide tight approx¬ 
imations for practical system dimensions of massive MIMO 
systems. Specifically, we focus on the following aspects: 

• We first characterize a lower bound on the feasible SINR 
targets for the ROBF algorithm. 

• We prove that the uplink and the downlink SINR achieved 
by the ROBF algorithm asymptotically converge to their 
target value in the large system regime. 

• We prove that the uplink and the downlink power allo¬ 
cations yielded by the ROBF algorithm asymptotically 
converge to the respective values of the CBF algorithm, 
hence making the ROBF algorithm optimal in the large 
system regime. 

We first start with the characterization of feasible SINR targets 
for the ROBF algorithm. 

A. Feasible SINR targets for the ROBF algorithm 

Throughout this subsection, we use the notations 7 to denote 
the vector of target SINRs as follows: 

7i = 7= [7i> • • • )7iv]- 

Similarly, we define the vectors 

— [^2,1) • • • 1 ^ — [^1? • • • 1 ^iv]? 

Mi = [Mu , ■ ■ ■, M= [Mi> • • • i Mat]- 

We now define the notion of feasible SINR target for the ROBF 
algorithm. Feasible SINR target implies that the following 
three conditions must be satisfied : 

[Cl] The iterations of the fixed point equation © must 
converge to a finite value //. This implies that given a 
target SINR vector 7 , there exists fi < 00 that satisfies 
the fixed point equation 

7 i,j = Vi,j. (19) 

Further from the property of the fixed point equation, for 
a given 7 there exists a unique fi satisfying ( fl9| (if the 
target SINR vector is feasible). Let us define such a pair 
of vector by {7 

[C2] For every pair of vectors { 7 , //} that satisfy ( fl9| ), the 
matrix I — TA must be invertible. 

[C3] The elements of the vector (I —TA) _1 p must be positive 

(in order for J to be real). 
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We focus on the conditions [Cl] and [C2] and defer [C3] to 
the end of the section. We first state the main result of this 
subsection and later on provide the proof. 

Theorem 1. Every target SINR vector 7 whose elements 
satisfy the condition 


K 


-Y- 

N t Ar 1 


7 i,k 


N K 


fc=1 


+ 7 i,k 


1 

W t 1 , 

n= 1 fc=l 1 + 
n^i 


c(n) 




c ( n ) 


< 1. (20) 


'7n,fc 


are feasible for the ROBF algorithm, where a ma x(n) =. 

We remark that in our result, the condition <© is only a 
sufficient condition. The proof of this theorem involves several 
steps which will be illustrated in the rest of this subsection. 
We proceed as follows. 

First, we establish a relationship between the two conditions 

[Cl] and [C2], 

Lemma 2. For every pair of vectors { 7 . //} that satisfy 
the matrix I — T A is invertible. 

Proof: The proof can be found in Appendix C, part I. ■ 
Lemma [2] implies that the condition [C2] is automatically 
satisfied for every pair of vectors { 7 , //,} that satisfy the 
condition [Cl]. Following the result of Lemma[2] we consider 
characterizing only the set of target SINR vectors 7 that satisfy 
the condition [Cl]. 

The exact set of feasible target SINR vectors satisfying the 
condition [Cl] is difficult to be characterized in closed form. 
Therefore, we only establish a lower bound on this set. In 
order to do so, we consider a modified system in which the 
inter-cell interference path loss coefficients are replaced by 


^"max (n) = sup Vn 7 ^ i. 

k 


( 21 ) 


where i = 1, ...,1V. Therefore, in the modified system 

Let us consider the ROBF algorithm applied to both the 
original and the modified systems. The fixed point equation 
for the computation of the uplink power allocation for the 
original system must satisfy the following equations (we use 
the superscript ”org” to represent the original system): 


org - org 

7 ij = ^.. 7 , 


v*, j, 


K 


org 

ViYk^ik 


k =1 
N K 


org - org 


N t 


YY 


org 

Oj,n,kM n 


■ y -1 , org - org 

1 + 07n,fc/Vfc m i 


71 = 1 k 

n^i 


+ 1 . 


( 22 ) 


(23) 


Similarly, the fixed point equation for the modified system must 
satisfy the following equations (we use the superscript ’’mod” 


to represent the modified system): 


^ , mod - mod 

7 i,j — Fi,k m i 


1 


_ ..mod 

u i,i,kH'ik 


N t 1 + a i>i , fc /r”° d m” od 

1 er max{n)p\ 


N K 

71=1 k—1 

n^i 


mod 

k 


c(n)p\ 


mod ^ mod 


1 . 


(24) 


(25) 


In what follows, we consider the ROBF algorithm applied to 
the modified system and characterize the feasible SINR targets 
corresponding to this system. The set of feasible target SINR 
of the ROBF algorithm applied to the modified system will act 
as a lower bound on the set of feasible target SINR of the 
ROBF algorithm applied to the original system. Intuitively, 
this is not hard to see. In the modified system, the path losses 
corresponding to the inter-cell interference links are scaled up 
to cr max (n). Therefore, the modified system represents a more 
interference limited regime as compared to the original system. 
Hence, any SINR feasible for ROBF applied to the modified 
system should be feasible for the ROBF algorithm applied to 
the original system as well. We will later on make rigorous 
arguments to prove that the above statement in Lemma [3] 

Proposition 2. The set of feasible SINR targets for the ROBF 
algorithm applied to the modified system must satisfy for i = 



Proof: The proof is provided in Appendix C, part II. ■ 
We now present the following Corollary. 

Corollary 1: 

In the large system regime, the set of feasible SINR targets 
for the ROBF algorithm applied to the modified system must 
satisfy for i = 1 ,N, 



+ 


N K 


N tYYi + - 


d n ) 


^n^k 


n— 1 k —1 
n^i 


“T n,k 


< 1 . 

(27) 


The arguments for the Corollary is also provided in Appendix 
C, part II. 

Finally, we show that the feasibility condition of ( |26| will 
act as a lower bound on the set of feasibile SINR targets for 
the ROBF algorithm applied to the original system. 

Lemma 3. Any target SINR feasible for the ROBF algorithm 
applied to the modified system is feasible for the ROBF 
algorithm applied to the original system. Hence the feasible 
SINR targets for the modified system will act as a lower bound 
for the feasible SINR target of the original system. 

Proof: The proof is provided in Appendix C, part III. ■ 
Lastly, we establish the condition [C3], This is rather 
straight forward and can be seen by the following steps: 























From Appendix C, we have established that for target SINR 
values that satisfy the condition in Proposition [2] p(T A) < 1. 
Therefore, using series expansion for the matrix (I — TA ) -1 
ED, we have 


(I-rA)-V= (]T(rA)')p. (28) 

j=i 

Since the elements of TA and p are positive, their sum will 
also be positive. Therefore, the condition [C3] holds true. 

We end this section by providing the feasibility conditions 
in two special cases in which the feasibility conditions are 
both necessary and sufficient: 

Example 1: Single Cell Case 

For the isolated single cell case, ([27]) reduces to the following 




Fig. 1. Downlink power Vs target SINR for K = 50 UTs per cell. 


B. Convergence of the Uplink and Downlink SINR 


When all the UTs are demanding the same SINR 7 *, = 7 Vfc, 
the set of feasible SINR is given by 7 < K y — lj 

The result can be interpreted as follows. In the case of a single 
cell, as long as N t > K , any finite SINR target is supportable 
by the ROBF algorithm (for the case when all 7 ?:,j are equal). 
In other words, when N t > K , the BS has enough degrees 
of freedom to sever all the UTs in the system. Note that ( |29| ) 
matches the feasibility conditions derived in lf20l for the case 
of a single cell system. 

Example 2: 2-Cell Wyner Model 

Consider a perfectly symmetric multi-cell system in which the 
path loss for the intra-cell links are equal to 1, and the path 
loss of the inter-cell links are equal to e. Every UT demands 
the same SINR target given by 7 . In this case, © reduces 
to 


N t \1 + 7 1 + 67 / 


(30) 


Next we focus on the achieved SINR in the uplink and 
downlink for the ROBF algorithm. We first start with the 
analysis of the uplink of the ROBF algorithm. Before we 
proceed, we make the following observations. Recall that the 
achieved SINR in the uplink for the ROBF algorithm is given 
by 


A&(/*) = 


|v ff h- ■ I 2 


E n-n,k 
(n,k)^(i,j) Nt 


>UDL\irH h . 
N t I V iJ Ul 


mi r 


(31) 


Note that in our notation A'A* (//). we have explicitly men¬ 
tioned the achieved uplink SINR as a function of the pa¬ 
rameters of the ROBF algorithm p (slightly deviating from 
the notation for the uplink SINR introduced in Section 0- 
Since the uplink receive filter in the ROBF algorithm ( fl3| are 
minimum mean square error (MMSE) form, the expression for 
the achieved uplink SINR can be given in alternate form as 


KM = v 


1 N t 




(32) 


This condition allows us to examine the dependency of the 
feasible SINR target on various system parameters. In partic¬ 
ular, we examine the dependency of the feasible 7 on N t by 
a simple numerical example. We plot the downlink transmit 
power (obtained by running the ROBF algorithm) as a function 
of 7 for different values of N t in Figure [I] We consider 
K = 50 UTs per cell and e = 0.5 It can be seen that beyond 
a certain cut off value of the target SINR, the downlink power 
grows unbounded. This is precisely the value of the target 
SINR at which the ROBF algorithm becomes infeasible. It can 
be verified that the cut off value of 7 is the one for which the 
condition in ([30]) is not satisfied. Finally, note that higher the 
number of transmit antennas per BS, higher is the cut off value 
of the target SINR. Once again, this is due to the availability 
of greater number of spatial degrees of freedom. In this cases 
when Nt > 2 K, any finite target SINR is achievable (note that 
2K is the total number of UTs in two cells). The feasibility 
conditions for the special case of two cell Wyner model was 
derived in CD as well. 


where = E i,j JvfKijhpij and S f = S M 

j. Similarly since the uplink receive filter <[8]i in 
the CBF algorithm are MMSE, the achieved uplink SINR for 
the CBF algorithm is given by 




3 


Vi, j. (33) 


Also, since the CBF algorithm is optimal. 


Vi, j. 


(34) 


Also, recall that the downlink SINR for the ROBF algorithm 
is given by 


A#(M) = 




E 




r H 




E 


n^i,k 3 I 


■ N 0 
(35) 


We now provide the following convergence result: 
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Theorem 3. In the large system regime, the achieved uplink 
and downlink SINR for the ROBF algorithm converge almost 
surely to the target SINR Ji j. Mathematically stating 


a$(m) 

K L » 


N t ,K-y oo 


> 7 i,j 


N t ,K-¥ oo 


> 7 ij 


Vi, 3, 
Vi,j. 


(36) 

(37) 


Proof: The details of the proof for the convergence of the 
achieved uplink SINR |36| can be found in Appendix D. 

Next we proceed to the convergence proof for the downlink 
SINR ( fT7j ). The proof utilizes the following lemma on the 
convergence result for the downlink interference terms. 


Lemma 4 . The downlink interference term corresponding to 
the intra-cell and inter-cell interference converge in the large 
system regime to the following: 


J2 \Vn,k h n,i,j\ 2 ~ J2 


^n,n,k^ttn 

N t 


(38) 


The lemma is proved in Appendix E. We now proceed 
to the convergence of the downlink SINR. Using Lemma [4] 
and Lemma uu it can be concluded that the downlink SINR 
asymptotically converges to 


A$(m) 


N+.K—tc 




N t ^2{n,k)j£(i,j) ^n,kG ni ijG nt n t k'm' n + Np 


(39) 


It can be easily verified from ( fl4| that the right hand side 
of ( [39] ) is equal to the target SINR 7 ij, thus completing the 
proof. ■ 


C. Asymptotic optimality of the Uplink and Downlink Power 
Allocation 

We now focus on the optimality of the uplink and downlink 
power allocation of the ROBF algorithm in the large system 
regime. 

Consider the Lagrangian of the downlink minimization 
problem in its two forms as in equations 0 and 0 - Our 
proof proceeds by plugging in the solution obtained by the 
ROBF algorithm (v, p) into the Lagrangian and examining 
its properties in the large system regime. 


Lemma 5. The following results hold true for the Lagrangian 
in the large system limit: 


lim L(v,p) = lim 


N+.K—t 


i,3 


PijN 0 


lim L(v, p) = lim > 

N—m*, Nt,K—too — J Nf 


N t ,K^ 


(40) 

(41) 


Proof: Please refer to Appendix F. ■ 

We can draw the following inference from the result of 
Lemma [5] 


Corollary 2: The uplink and downlink power allocations 
yielded by the ROBF algorithm are equal in the asymptotic 
limit: 


lim 

N t ,K-too 




hm 

N t ,K—fOO 

1,3 


Pi,jN 0 

Nt 


(42) 


We know that the uplink power allocation of the CBF algo¬ 
rithm satisfies 

\ ^—LI —- = minmaxLfw, A) (43) 

N t w a v ' 

h3 


If we show that in the large system regime 


lim > 

N t ,K^t 00 ' 


Pi,jNp 

N t 


lim min max L (w, A), 

N t ,K—too w A 


(44) 


and the duality gap is zero, then this implies that the so¬ 
lution provided by the ROBF is optimal in the asymp¬ 
totic domain. In other words, the optimal downlink power 
(which is the solution of the primal problem) is equal to 
limjv ti if_).oo min w maxA L{ w, A). In order to do so, we prove 
the following result: 


Lemma 6. In the large system regime, the sum of uplink power 
allocation of the ROBF algorithm converges to the sum of 
uplink power allocation of the CBF algorithm. 


E Pi,} „ Aj 

Nt ~ 4^ Nt 


Vi,j. 


(45) 


Proof: The result is proved in Appendix G. ■ 

Consequently, the solution provided by the ROBF algorithm 
is optimal to limjvt,jr->.oo max w minA L(w, A). Also, from 
the result of Corollary 2, it follows that the downlink power 
allocation of the ROBF algorithm is also optimal in the large 
system regime. This concludes the proof. 


IV. Incorporating Individual BS Transmit Power 
Constraints 


In this section, we consider a practical constraint in MIMO 
multi-cell systems, namely the individual BS transmit power 
constraints. Recall that the basic optimization problem consid¬ 
ered in this work stated in ([ 2 ]) does not impose this constraint. 
In what follows, we propose a heuristic extension of the ROBF 
algorithm to incorporate this constraint. 

The optimization problem in Q along with the individual 
BS transmit power constraints can be stated as follows: 


min 

W i,j Mi,3 

S.t. 


J2 

i,3 

Vi, j, 

^ ^ < -^ 2 ,max 


Vi, 


(46) 


where Pi max denotes the peak power of BSj . In what follows, 
we first state the heuristic extension of the ROBF algorithm to 
incorporate the individual BS transmit power constraints, and 
then provide the main intuition behind the development of the 
algorithm. 
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Algorithm 3 (ROBF with individual BS transmit power 
constraints). Perform the following steps. 

1. Initialize r = 0, and cti(O) > 0, Vi. 

2. Starting from any initial p® j(r) > 0 Vi,j compute the 
uplink power allocation as Pij(r) = linij^oo p\ ■ (t), 
where 




7 i, : 




-A(t) 


Vi, j 


(47) 


and m\(r) is evaluated as m\[f) = limp-^ m*’ p (r) 
(initializing with any to*’°(t) > 0,Vi) 

= (48) 


where a, : are the Lagrange multipliers associated with the indi¬ 
vidual BS transmit power constraints, and a. = [a i,..., ajv] T 
(Ahas the same interpretation as in ([3])). From ([T}, we can 
rewrite 0 as 

L'{ w, A, a ) = L(w, A) + ^ Ctj [ ^ W U W i,j ~ P I max • 

* i 

(54) 

Using duality theory, the solution of the dual problem can be 
given by max a ^ min w Z/(w, A, a). Let us denote 

g(a) = max min 2/(w, A, a). 

A w 


3. 


Vi,n,kPn,k( T ) 


-1 


v JV * 7X 1 + Cr i,n,fc^,fe(' r )w- 
Set f/ie receive uplink filter as 
= 

1 / Pn,k( T ) 

7^ 


,p-l(r) + 1 + a i( r ) 


(49) 
-1 




n,k 


N f 


hi,n,fehj n J. + (1 + Cti(r))I^ 

(50) 


Our approach proceeds by showing that j w 7j w i,j ~ p i, max 
is a sub-gradient direction of the function g(a). Then, the 
dual problem max a [max^ min w L'(w, A, a)] can be solved 
by updating the Lagrange multiplier a in the direction of the 
sub-gradient Il24ll . 

In order to do so, consider two vectors a ' 11 and <y r2> . Note 
that 

g(a^) = maxminL'(w, A, A (1 \ a^), 

A w 


4. Set the transmit downlink beamforming vectors as 

Vj j(t) = \j ^i,ji where the scaling factors Sij(r) 
are calculated as in 0- 

5. . Set t = t + 1 and update afir + 1) as 


otfir + 1) = 



+ ^(E v 4i( r ) v *.i( T ) 

j 



Vi, 

(51) 


where £ > 0 is a small step size. If cti = 

0 ,Vi € {1, then terminate. Else if for all 

i £ {1,..., TV}, for which cti > 0 if 


E v 5( r ) v L/( r ) - p ,max 
3 


< 6 , 


(52) 


where we use (w^ 1 ), A* 1 -*) to denote the solution to 
max\ min w L'( w, A, a( 1 )). Similarly, let (w^ 2 ), A (2 ^) denote 
the solution to max^ min w L'( w, A, a£ 2 )). Consider g(u W). 
Firstly, it can be noted that 

min L'(w, A, a^) < L'(mv( 2 \ A, a^), (55) 

W 

and hence 

maxminL'(w, A, a*- 1 -*) < maxL 7 (w (2 \ A, (56) 

A w A 

Using the definition of ^(aA 1 )) in ( |56| . we conclude that 

3(q 4)) < maxL , (w* 2 *, AjQ* 1 -*). (57) 

A 

Using 0 in the right hand side of ( |57] , we obtain, 

g(a (1) ) < 


where S > 0 is a small non zero quantity, then terminate. 
Else return to Step 2. 


We now proceed to provide the main intuition behind this 
algorithm. 

A. ROBF with individual BS transmit power constraints: In¬ 
tuition 


maxL(w (2) , A,a (1) ) + ^“! 1) [E( w l) 


H w?} - p , 


(58) 


Adding and subtracting the term <:7 P> ~ 

Vi,max to the right hand side of ( |58| , and rearranging the 
terms, we obtain 


,(2)\J?,„(2) 


Consider the optimization problem in |46|. In order to 
solve this problem, we proceed by considering the Lagrangian 
associated with ([46]), given by 


L'(w, A, a) = 0 wf.'.vv,., 
ij 


E 


Nt L 7 i,j 


^ ^ |w njfe h n)iJ | iVo 

(n,k)^(i,j) 

“t“ ^ ' OH ^ W,; .y Pi,max , 

* 3 


( 53 ) 


g(ot P 7 < maxi'fw^ 2 *, A,a^) 


E( 


„(!) 


a; — a 


( 2 ) 


[5>$> 


"w! 2 ? - p, 


i,3 


v(« ( 2 ) )+eK ) -^ ( 2 ) )[E( w S) 


- R 


(59) 


where in ([59], we have used the definition of g(oc^). Conse¬ 
quently from ( |59| , it can be concluded that Sj( w i^)’ H ' w i j* 
P.max is a sub-gradient direction of the function g(a) (by the 
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definition of a sub-gradient). Therefore, updating the Lagrange 
multiplier cti in the direction of the sub-gradient solves the 
dual problem max Q [max^ min w L'(w, A, a)] J24|. There¬ 
fore, a,; must be updated as 


adT + 1) = 


a *( T ) +C(5Z W 5( T ) W id( T ) - P I max) j Vf, 


(60) 


where ( > 0 is a small step size. 

Further, for each value of a(r ) one needs to solve the 
problem max^ min w Z/(w, A, a). This can be solved by 
repeating the steps of the CBF algorithm for non-zero values 
of a. 

Recall that all the above arguments were made considering 
the case of finite system dimensions. We now follow the 
same approach as in the development of ROBF algorithm, 
i.e., in the large system domain, utilize RMT results to obtain 
asymptotic approximations of the quantities involved. The 
ROBF algorithm with per the BS peak power constraints 
developed in this section is obtained by following this idea. 
The theoretical analysis of the performance of this algorithm 
along with the proof of convergence will be a topic of future 
research. Herein, we resort to the numerical results stated in 


Section VI in order to show the performance as well as the 
convergence. 


V. Impact of Imperfect CSI and Pilot 
Contamination on the performance of ROBF 
Algorithm 

In this section, we investigate the impact of CSI estimation 
errors and pilot contamination on the performance of the 
ROBF algorithm. Throughout this section, we assume that 
the slow fading co-efficient (path loss information) can be 
accurately estimated at the BS (since they remain constant 
for a long period of time, they are easy to estimate, see 
for e.g. l25l ). Further, for the fast fading co-efficients, we 
assume reciprocity between uplink and downlink channels, and 
consider the time division duplexing (TDD) model of channel 
estimation, i.e., estimation via uplink pilots. 

A comprehensive design of the optimal beamforming vec¬ 
tors in the presence of CSI estimation errors and pilot con¬ 
tamination issue is out of the scope of this paper. Alternately, 
we take the following approach: First, we assume that the 
BS treats the CSI estimate as the true CSI and implements 
the ROBF algorithm directly. In this case, we derive the 
asymptotic equivalent of the downlink SINR achieved by the 
ROBF algorithm. Using numerical results, we investigate the 
impact of CSI estimation errors and pilot contamination on 
the performance of the ROBF algorithm both in terms of 
the achieved SINR and downlink power. Then, exploiting 
the fact that the BS has the accurate knowledge of the 
slow fading co-efficient, we propose a heuristic adaptation of 
the ROBF algorithm in the presence of imperfect CSI and 
pilot contamination named as the modified ROBF (MROBF) 
algorithm. We show that under the massive MIMO regime, an 
algorithm in which parameters can be computed based on the 


channel statistics (rather than the fast fading CSI) such as the 
MROBF algorithm is more robust to CSI estimation and pilot 
contamination effects. 


A. CSI Estimation 

We now describe the CSI training phase for the estimation 
of the fast fading co-efficients. Let T be the length of the 
channel coherence interval, a part of which is dedicated for 
CSI estimation, and 7), be the number of symbols used for 
pilots. Therefore, the UTs in every cell i transmit Tj r mutually 
orthogonal pilot symbols to their respective BSs during the 
training phase. We represent the pilot sequences used by the K 
UTs in each cell by the matrix \j € C TjrXK , (Tx r > K). 
The pilot sequences are repeated in each cell, hence leading 
to the issue of pilot contamination. The matrix $ satisfies 
$ h $ = I k . 

The signal received during the channel training phase de¬ 
noted by Y i Tr £ C ' Vl x K can be written as 


Y i>Tr = H m $ t + ( H i,„) $ T + N, (61) 

n^i 

where H ijn = [h i>Tli i,..., h itri:K ] and N £ C' ViXTt ' with 
i.i.d. CAf(0 ,1) elements represents the noise during channel 
training phase. The MMSE estimate of h, „ given Y, Tr can 
be given by 1261 


N 


hi,n,k 


(62) 

where = a i,n,k (Ef=l + i^) 


. It can be veri¬ 


fied that h, „ k is distributed as 


N 


-IN 


hj : ra,fc ~ CM I 0, n I ^ „ 


(63) 


s.6=l 


For notational simplicity, let us denote cr in / i: = 

°ln,k (E£Li <Ti,b,k + 3 ^) ■ Note that 


a i,n,k ~ a i t n,k a i,n,k- 


(64) 


ROBF Algorithm With CSI Estimates 

Throughout this subsection, we assume that the BSs assume 
the CSI estimates to the true channel vector, and implement 
the ROBF algorithm. Note however that the computation of 
fiij and Sij in eu and < |T~4| > do not require the estimates of 
the fast fading co-efficient, and hence they can be implemented 
directly. The uplink receive filter with imperfect CSI denoted 
by vfj can be formulated as 
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The dowlink beamforming vector denoted by vfj can be then 
computed as 



We now investigate the achieved SINR in the downlink 
under the ROBF algorithm with imperfect CSI. The expression 
for the dowlink SINR can be given as in ( |35| , by replacing 
Vj j with . The asymptotic equivalent of the dowlink SINR 
can be obtained by analyzing it in the large system regime. 


Theorem 4 . The achieved SINR in the downlink by the ROBF 
algorithm in the presence of imperfect CSI converges almost 
surely to the right hand side of <HD in the large system regime, 
where the term G n ^ n ^ is defined as 


1 v - '' N 

A V j * 1 &i,n,kBn,k 

C n ,n,k = , , >. - 2 - and t, it k = -. 

1 “f C,n,kdn,n,ktri n &i,i,k 


( 68 ) 


Further, raff can be computed as the solution to the fixed 
point equation 



£,n,kdn ,n,k 

1 T fn,k.dri,n,ktit e f 



(69) 


BS cannot distinguish between the channels of UT, j and 
UT n j due to the pilot contamination effect. This implies that 
naive application of the ROBF algorithm may significantly 
underestimate the interference arising out of the UTs that reuse 
the same pilots. The algorithm performance can be enhanced 
by carefully accounting for these issues. This is indeed the 
main intuition behind the MROBF algorithm. 

In the MROBF algorithm, we consider that the BS does 
not alter the structure of the beamforming vector, i.e., the BS 
retains the RZF beamforming vector as in ( |67] i. Nevertheless, 
the performance gain can still be obtained by redesigning the 
uplink power allocation, and computation of S. We redesign 
the uplink power allocation on similar lines as that of the 
ROBF algorithm. 

First note that in the case of perfect CSI, the achieved 
uplink SINR in the asymptotic limit is given by <7i } i,jPi,jrhi. 
Further, the uplink power allocation is chosen to satisfy 

7 i,j = or 


Bid — 


7 i,j 


Vi, j. 


(77) 


We use a similar argument for the computation of the uplink 
power allocation in the case of imperfect CSI. In order to 
do so, consider the uplink SINR with uplink receive filter 
formulated as in 


and ( m' n ) es ' can be computed from rrfff as 

Kf) 2 


«r = 


Nt ^k=l (X+v n ,n,kin,k™“ 


l _ j- 

1 N, 2^k= 


Proof: The proof is provided in Appendix H. 


(70) 


B. Modified ROBF (MROBF) Algorithm 

In this subsection, we propose a heuristic adaptation of 
the ROBF algorithm addressed as Modified ROBF (MROBF) 
algorithm which is designed to accommodate the effects of 
imperfect CSI and pilot contamination. 

First, note that naive application of the ROBF algorithm 
may not yield good performance in the presence of imperfect 
CSI and pilot contamination. The following reasoning provides 
an intuitive understanding for this: Consider the interference 
arising from the signal of UT n (n i,k ^ j) at the UT,;^ 
(which does not use the same pilot as UT;.j), i.e., the term 
|w„ ifc h niiJ | 2 . Following the derivation of Appendix H, it can 
be verified that this interference terms converges to 


w. 


. 1 2 — ^ n,k ~ n -2 / - / \est p 

n,k a n,i,j\ ,r <-'n,rt,fc ( - r n j r li fc \ m n ) £>n,k- 


N< 


(75) 


Now consider the interference arising from the signal of 
UT n j, (n i) at the UT, ? (which reuses the same pilot 
as UTij), i.e., the term |w n jh n i j| 2 . It can be verified that 
asymptotically, this interference terms converges to 

&n,j ( a n,i,jO' n n jG ni n,jtn n • (76) 


We note that the two interference terms are significantly 
different. In particular, CD is scaled by a factor of . 
(and for large N t , this has a very low value). However, the 
term in © is not scaled. This is due to the fact that the 


(A£( M ))e* = 


FLd. |v est h- ■ I 2 


(n,k)ji(i,j) N t 




hJ 1 


(78) 


Let us examine the uplink SINR in the large system domain. 
Using the derivation similar to Appendix H, by replacing the 
individual terms of ( |78| by their asymptotic equivalents, the 
uplink SINR in the large system limit can be approximated 
by (|73lP1 For convenience, let us denote that denominator of 
f73j ) by X^y mp . Following the same approach as in ( |77j ), we 
compute the uplink power allocation as the solution to the 
following set of equations: 


Bid — 


Hd 


asymp 


{vi,id G i,id m i St Y 


Vi, j- 


(79) 


The pij that satisfies the set of equations ( |79[ i can be com¬ 
puted using the following iterative method: Starting from any 
initial pf 1 > 0 Vi, j the uplink power allocation is given by 

Bid = ^ m t-kooBi,j> where 


u t+1 = 

B%,j 


N- (J UL ) 4 
h,3 V-^asymp ) 


' T i,id G i,id m '- 


(mf 1 ) 4 ) 2 


Vi, j- 


(80) 


In m , (T^ymp) 4 , GC :) and (mf 1 ) 4 denote the respective 
quantities computed at p\ j ; Vi, j. We observe that numerically 
that the iterations of ([80| converges to the solution of ( |79| ) (see 
numerical results in Section |VI| ). 

Finally, 5 can be computed as a solution to the following 
linear equations 


S = A -1 7 


(81) 


3 Note the this approach does not constitute a formal proof of convergence 
of the uplink SINR. Since the algorithm developed itself is heuristic, the proof 
of convergence has been omitted. 
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(A£(M))est 


where 




12n=lS n ,j {^ n ,i,j^' n n jG nyn jm^\ + ^ En=l E*=l 3n,kd'n,n,kGn y n,k( m n) e> ' t Bn,k 

n^-i \ J 


ddn,k ) Gn,n,j‘kkln ■ 

En=i MnJ En=l Efc=l Ibi.k&i.ij&i,i ,j (K ) est B; fc + a^jG, 2 . -(m^ 

where 

^n,/c = &i,n,k T ^ {G,k&i,n,k& i y i y kGi y i,k^i )~ ‘^ l L* J i,k(.&i,n,k&i y i } k) Gi y i y klTl^ . 


(71) 

(72) 

(73) 

(74) 



where the matrix A £ C m " xM is defined as 



where each submatrix A ,J £ C KxK is given by 


(82) 


' (vi,i,jGi,i,jrnf t ) 2 , 


= 


— 1 a. /^t2 / - / \est d 

JSf t u n,n,k^ r n ^ n ^k\ ,,l 'n) D n,lc? 


) 2 . 


n = i, k = j 

n^i,k = j 

n = i,k^j 

and n ^ i,k / j. 

(83) 


VI. Numerical Results 

In this section, we present some numerical results to demon¬ 
strate the performance of the ROBF algorithm in a massive 
MIMO setting. 

We consider a hexagonal cellular system consisting of 2 
cells as shown in Figure [2] We use the distance dependent 
path loss model in which the path loss from BS i to UTj is 
given by 

do 


where d ly jj,- represents the distance between BS i to UT^fc. 3 
represents the path loss exponent which is taken to be 3.6 in all 


the simulation scenarios, do represents the channel attenuation 
at reference point and is taken to 10 -3 53 . Location of the UTs 
are obtained by generating uniform random numbers inside 
each hexagonal cell. The distance d m ; n < d l:j k < d max , 
where d m i n and d max represent the minimum and the max¬ 
imum distance between the UTs to the BSs of their respective 
cells. In our simulation, d m i n = 20 m and d max = {500,1000} 
m depending on the distance between the two BSs. The noise 
power is taken to be —104 dBm over the operating bandwidth. 
All numerical results are plotted by varying the positions of 
the UTs inside the cell over 500 iterations. 


First, we examine the performance of the ROBF algorithm 
in satisfying the UT SINR constraints. Accordingly, we plot 
the variation of the uplink and downlink SINR (averaged 
across the UTs) as a function of the number of antennas per 
BS for 1000 channel realizations in Figure [3] Herein, K = 50 
UTs per cell and target SINR is 3 dB per UT. The horizontal 
line represents the target SINR which is 3 dB. The bubbles 
represent the average achieved SINR values (averaged over the 
channel realizations), and the vertical lines around this bubble 
represent the variation of the achieved downlink SINR around 
the average value. It can be observed that even for moderate 
number of antennas, e.g. 60 antennas per BS (comparable to 
the number of UTs), the target SINR constraints are satisfied 
for almost every channel realization (since the fluctuations are 
small). This implies that the ROBF algorithm nearly optimal 
under this setting in terms of satisfying the SINR constraints. 


Next, we investigate the downlink power of the ROBF 
algorithm and compare it wit h the CBF algorithm and ZF 
beamformiiH' 


PfWG ,, where w?} is unit 

L iJ t iJ L iJ 


(denoted by 

norm vector Jr 1 Note that our simulation setting, the BS per¬ 
forms ZF beamforming to null the interference UTs in both the 
cells (and not merely the UTs in its own cells). After nulling 
the interference, BS performs appropriate power allocation in 


4 Note that for eigen beamforming and RZF, after fixing the beamforming 
direction, the power allocation needed to satisfy the SINR constraints does 
not yield to a simple structure (as in the case of ZF). In order to avoid the 
additional complexities associated with power allocation, only ZF is used 
for comparison purposes. Nevertheless, we expect the gains of the ROBF 
algorithm to hold true with respect to the other beamforming techniques as 
well. 
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Fig. 3. Fluctuations of the downlink SINR (ROBF algorithm) around theFig. 4. Comparison of downlink power per UT as a function of the number 
target value. K = 25 UTs per cell and target SINR = 3 dB per UT. of antennas per BS. K = 25 UTs per cell and target rate = 3 bits/s/Hz 

(log(l + 7t,j)) P er UT. 




15 

No. of UTs per cell 


Fig. 5. Comparison of downlink power per UT as a function of the number Fig. 6 . Comparison of downlink power per UT as a function of the K. 
of UTs per cell. Nt = 60 antennas/ BS and target rate = 3 bits/s/Hz Settings identical to Figure [5] 

(log(l + 7 i j)) per UT. D represents the distance between the two BSs. 


order to meet the target SINR constraint of the UTs 

pZF _ Ji,jNp 

iJ 

We plot the downlink power per UT (i.e. Sumdo—p ? w 5) afj 
function of the number of antennas per BS in Figure [4] Herein, 
K = 25 UTs per cell and target rate = 3 bits/s/Hz (log(l + 
7 ij)) per UT. The following conclusions can be drawn. First, 
the downlink power expended by ROBF very closely matches 
that of CBF, indicating that the ROBF is nearly optimal in 
terms of minimizing the downlink power as well. Next, it can 
be seen that for moderate number of antennas per BS, e.g. 
50 — 100 antennas, ROBF provides substantial power gains as 
compared to ZF beamforming. This result highlights the gains 
obtained by optimizing the power allocation (as compared to 
naively nulling out interference to all the UTs in the system). 
When Nt is very large, the scaling of ^ starts to play a 


dominant role, and ROBF ceases to provide substantial gains 
over ZF. 

Finally, we examine the efficacy of ROBF algorithm in 
terms of its ability to support greater number of UTs per cell. 
This is accomplished by fixing the number of antennas per BS 
to 60, and plotting the downlink power as a function of the 
number of UTs per cell in Figure [5] It can be seen that beyond 
a certain number of UTs, the downlink power corresponding to 
the ZF beamforming becomes unbounded. In fact, this happens 
at K = 30 UTs per cell (note that N t = 2 K at this point 
indicating that the BS has used up all its degrees of freedom). 
At the same time, the downlink power of the ROBF algorithm 
grows unbounded at a much later stage, i.e. around 55 UTs per 
cell. Moreover, we can also see that as D (distance between 
the BSs) increases from 500 m to 1000 m, the transmit power 
decreases as expected. This is due to the reduced effect of 
inter-cell interference. Lastly, to quantify the transmit power 
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gains obtained by the ROBF algorithm over ZF, we zoom into 
the previous plot in Figure [6] Once again, as noted before, 
ROBF algorithm provides substantial power gains compared 
to ZF, a gain of 5 dBm per UT for K = 25 UTs per cell. 


ROBF with Individual BS Transmit Power Constraints 


In this subsection, we illustrate the performance of the 
ROBF algorithm with the individual BS transmit power con¬ 


straints developed in Section IV In particular, we show the 
convergence as well as the performance results of this algo¬ 
rithm. For illustration purposes, we consider a 2 cell scenario 
and an identical system set up as in the previous case. The 
peak power of BSi is set to 20 dB and the peak power of 
BS 2 is set to 10 dB. We run the algorithm of Section IV 
until the stopping criteria mentioned in Step 5 is satisfied. 
The transmit power per BS is plotted against the number of 
iterations in Figure [7] The horizontal dashed lines in this figure 
correspond to the transmit powers obtained by running the 
original ROBF algorithm (without the individual BS transmit 
power constraints). It can be seen that for this case, the peak 
transmit power of BS 2 is greater than 10 dB, and hence its 
peak power of constraint is violated. The solid lines represent 
the transmit powers obtained by running the ROBF algorithm 
with the individual BS transmit power constraints. It can be 
seen that the transmit powers of both the BSs converge within 
a few iterations. Further, the transmit power of BSi is equal 
to 10 dB, and hence its peak transmit power constraint is 
respected. The transmit power of BS 2 is greater for this case 
as compared to the original ROBF algorithm. Further, we 
observed from the numerical results that the total downlink 
power of the system (summing the transmit power of both the 
BSs) obtained by running the original ROBF algorithm is 14 
dB where as its value for the case of ROBF algorithm with 
the individual BS transmit power constraints is 14.65 dB. This 
illustrates the fact that the minimum downlink power with the 
individual BS transmit power constraints is greater than the 
minimum downlink power of the original optimization in |2|. 


Impact of Pilot Contamination Effect 

We finally investigate the impact of pilot contamination 
effect on the performance of the ROBF algorithm, and the 
improvements obtained by the MROBF algorithm. We con¬ 
sider an identical system set up as in the case of perfect CSI. 
Since our main objective is to characterize the performance 
loss due to the pilot contamination effect, we ignore the errors 
associated with the CSI estimation process, i.e., Pj T in ( |62| ) is 
set to a very high value. 

First, we compare the achieved SINR in the downlink by 
the MROBF algorithm and the CBF algorithm in Figure [8] 
For the implementation of CBF, we consider that Algorithm 
1 is implemented with all the CSI values replaced by their 
estimates (i.e. ii l n .k replaced by h Wi,n,k computed as 
in (|62|)). It can be seen that the MROBF algorithm very closely 
matches the target SINR values, where as the CBF algorithm 
does not achieve the target SINR. This is due to the fact that in 
the MROBF algorithm, parameters can be computed based on 



Number of Iterations (t) 


Fig. 7. Iterations of the downlink power per BS. Nt = 100 antennas/ BS, 
K = 50 UTs per cell, Pi, ma x = 20 dB and P 2 ,max = 10. Target rate per 
user = 3 bits/s/Hz. 


the slow fading co-efficient, which can be estimated accurately 
at the BS. However, the CBF algorithm requires knowledge of 
the fast fading co-efficient, whose estimation suffers from the 
pilot contamination effect (in addition to the CSI estimation 
error that has been ignored in our numerical results). 

Finally, we provide a comparison of the downlink power 
reduction provided by the MROBF algorithm. Notice that the 
CBF and ROBF do not achieve the target SINR. For the sake 
of comparison, we devise an algorithm that achieves the target 
SINR. One can think of retaining the uplink power allocation 
of the ROBF algorithm, and only adapting the computation 
of 5 as in such that the target SINR is achieved in the 
downlink. We address this algorithm by the name ROBF with 
downlink adaptation. We plot the transmit power per UT as a 
function of the target SINR for both the MROBF algorithm 
and the ROBF with downlink adaptation in Figure [9] It can 
be seen that the transmit power per UT is significantly lower 
for the MROBF algorithm. Moreover, beyond a certain target 
SINR value, the downlink power for the ROBF with downlink 
adaptation becomes very high, implying that the target SINR 
cannot be supported by this algorithm. In contrast, the MROBF 
algorithm can support a greater range of target SINR values. 
This is due to the fact that the ROBF with downlink adaptation 
algorithm significantly under estimates the interference arising 
out of the UTs which reuse the same pilot symbols, where as 
the MROBF algorithm accurately accounts for this quantity. 

VII. Conclusion 

In this work, we formulated a decentralized multi-cell 
beamforming algorithm for massive MIMO systems with the 
objective of minimizing the aggregate transmit power across 
all the BSs subject to UT SINR constrains. The algorithm 
requires only locally available CSI at the BSs, and some 
statistical side information of the channel gains to other UTs, 
and incurs a lower burden both in terms of the information ex¬ 
change between BSs and the number of computations. Further 
we proved theoretically that this algorithm is asymptotically 
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Fig. 8. Achieved SINR Vs Target SINR for MROBF and CBF algorithmFig. 9. Comparison of downlink power per UT Vs Target SINR for the 
in the presence of pilot contamination effect for Nt = 60 antennas/BS andMROBF and ROBF algorithm with downlink adaptation. 

K = 30 UTs per cell. 


optimal for a large number of BS antennas and UTs. We also 
characterized a lower bound on the range of SINR values 
for which this algorithm is feasible. We confirmed using 
numerical results that the algorithm is nearly optimal in the 
regime of massive MIMO systems. It also provides substantial 
power savings as compared to zero-forcing beamforming when 
the number of antennas per BS is of the same orders of 
magnitude as the number of UTs per cell. We proposed a 
heuristic extension of this algorithm to incorporate a practical 
constraint in multi-cell systems, i.e. the individual BS transmit 
power constraints. We verified the convergence as well as 
the performance of this algorithm numerically. Finally, we 
investigated the impact of CSI estimation errors and pilot 
contamination on the performance of the ROBF algorithm. 
Further, we proposed a heuristic modification of the ROBF 
algorithm, and showed with the help of numerical results 
that such an algorithm is more robust to the effects of CSI 
estimation errors and pilot contamination (as compared to an 
algorithm that uses fast fading CSI values to compute the 
system parameters). 


Appendix A: Some Relevant Results 

Lemma 7. (?[27l?) For two matrices X and Y. if() < X < Y. 
then p(X) < p(Y). 


Lemma 8. (Equation 2.2, l\28j ) Let A. be a Hermitian 
invertible matrix of size N X N, then for any vector x G C N 
and scalar t £ C for which A + txx h is invertible, 


(A + rxx ff ) 1 = A 1 


A _ 1 rxx i? A _1 
1 + rx H A _1 x ’ 


and 


x H (A + rxx fl ) 1 


x H A~ 1 
1 + tx h A -1 x 


Lemma 9. (Lemma 2.6, /29l/ > Let x, y ~ CAf(0, 77 Tv) G 
C N , A G (£ NxN Hermitian matrix, and that x and y are 


mutually independent, and independent of A. Consider in > 
2. Then there exists a constant C m independent of N and A 
such that 


E 


x Ax —tr( A) 


< 


C„ 


Nm/2 I 


This implies by the Markov inequality and the Borel Cantelli 
lemma for ||A|| < 00 , that 

x H Ax - -^tr(A) - as '- > 0. 

N V ' N^r OO 

Additionally, 

H A a - s - n 

x Ay- > 0. 

N—foo 


Lemma 10. (Lemma 2.6, ( 1251? ) Let z G C + with v = Im(z) 
and A and B are N x N matrices with B being Hermitian, 
rGl, and q G C N , then 

|fr((B — zT)- 1 - (B + rqq ff -^I)" 1 )A)| < M. 

V 

Theorem 5. ( M8\D Consider the matrix B y, = XAr t Tjv t X^ t , 
where Xjv t = -j=Y Nt G (£ NtXNh w ith entries Y Nt (p,q) ~ 
Cff(0, 1), and the matrix Ttv, a non random diagonal matrix 
given by TV t = diag(t\ ,..., Ink) G ^ NKxNK _ p e f m (z) = 
-h-tr( B^ + zl) _1 , z > 0, Then, m(z) — m(z) -—-» 0, 

4 Af t ,A'->oo 

where fh(z) can be evaluated as the unique solution to the 
fixed point equation 

( 1 ^5 ti \ 

Further from m, under the same assumptions as Theorem 
|5j let m'(z) = j^tr(Bjv t + zl)~ 2 . Then, 

m! (z) — m(z) -—-> 0 (84) 

Nt ,K—f 00 
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where m'(z) can be evaluated as 


m! (z) 


m 2 (z) 

m 2 (z) ^NK if 

1 N t Z^i= 1 |14.tifn(4P 


(85) 


At this point, we remark that all the above results also hold 
for a more general non-Gaussian vectors/ matrices satisfying 
some moment conditions. 


Lemma 11. (Lemma 1 3301 ) Let a n, Ay, x n and yjy denote 
four infinite sequences of complex random variables indexed 
by TV. If cia i x b n and xn ~ UN and if | a at |, | J/1 tv 1 
and/or |ajv|, y | y” 1 are uniformly bounded above over TV, then 
cln/xn bisr/yN- The requirement of uniform bound can be 
relaxed to boundedness almost surely. 


Appendix B: Proof of the Fixed Point Equation 

In order to prove the convergence of the iterative equation 
we use the arguments of standard function ED- a 
A'— variate function g(x) = [g 1 (x),... ,g^(A] £ for 
x £ C A is said to be standard if it fulfills the following 
conditions: 

• Positivity: g(x) > 0, for x > 0. 

• Scalability: For (3 > 1, /3g(x) > g (fix). 

• Monotonicity: For x' > x, g(x') > g(x). 

Now consider the following iterative algorithm given by, 

x t+1 = g(x 4 ), t > 1 (86) 


If the A'—variate function is standard, then it ensures the con¬ 
vergence of [86) to its unique fixed point solution x = g(x), 
if the solution exists. 

Consider the iterative equations in eu which can be 
represented as 

= fidtf) (87) 

where p\ = [p\^. .., p\ K } T and p { = [p\,..., /4J r and 

kAh) = 1 ( 88 ) 

Let us also define the TV AT— variate function 

f(/x‘) = [f!(p‘).f w (/if eR™. (89) 


The existence of the fixed point to the equation p = i(p) 
follows from Lemma 1. We now prove that the TV A'— variate 
function f (p) is a standard function and hence the iterations 
of ( [87] ) converges to its unique fixed point. 

In the subsequent part of this proof, we introduce the 
notation fh\(z) to denote the solution of the fixed point 
equation 


to 4 0) 


V 1 n—1k=1 





With this notation, the solution to the fixed point equation of 
0 can be written as m 4 ( 1). For notational convenience, we 
also drop the superscript t. 

Positivity: The positivity result follows directly since rhi( 1) 


is positive whenever p> 0. Hence fij(p) > 0. 

Scalability: Let us consider the difference between the 

following quantities. 


PfipA) - fid Ah = 


hid 


7 id 


7 id ( P 1 \ 


' 1 , 1,3 \ m. 


(i) 


(i) 




(91) 


where (3 > 1 and m[ 1 \ 1) and m( 2 ^(l) are the unique solutions 

to ( fl2l > evaluated at p and (3p. In order to evaluate frh (1), 

( 2 ") 

we go back to the definition fh\ ' (1) evaluated at (3p. 


rn?\ 1) = ^tr(H,/3MHf + I)" 1 x Jm, (£) (92) 


where M = diag[/x]. Clearly 1//3 < 1. m| 2 (( 1) can be 
evaluated as the solution to the fixed point equation © 
evaluated at the point 1/(3 and then scaling the result by (3. 
From [9l) and [92) , it can be concluded that in order to prove 
the scalability result, it is sufficient to show that mfiz) is a 
decreasing function of z. 

In order to prove the same, let us consider an extended 
version of the channel matrix which is constructed as follows. 
Defining, R iJifc = a i , J yl Nl £ R N±xNt and Rfu fc = 
f Xi.j,k~3-N t .L £ R N±LxNtL . The matrix Hf is constructed as 
follows: 


1 


Hf, = y- 
J s/L L 


id. 

Hf = £ 


-pl/2 t> 1/2 


N+LxNKL 


x N t LxKL 


(93) 


where the matrix G C NtxK , whose elements 

X-^j(p,q) ~ CAf(0, ^). Also, let us define the following, 



j ^, Xij 


\ L \ 1 

A i,l> ’ ’ * > A i 


1 T 


K 


\)Lx 1 


t)KLxl 


diag ( 


A 1 ? * • • 5 A jV 


kNKLxNKL 


(94) 


and let Qi(z) = (HjAHf + zI^l) 1 and Q,f(z) = 
(Hf A L (H L )H + zlYti) -1 - Let us denote 


™f(z) 


1 

NtL 


tr(QfW). 


(95) 


It can be verified that for any fixed TV, TV t and I\, the following 
limit holds, 


m i( z ) ~ m i{ z ) - as - > 0. (96) 

L—foo 


Now consider the difference between the following two quan¬ 
tities, for any z^ > z\ > 0 and for any positive L, we have 

mf (zi) - mf(z 2 ) 

= - QfM 

(a) 

> 0 (97) 
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where inequality (a) follows by the following identity: For 
invertible matrices A and B, 

A " 1 -B " 1 = -A-^A-BjB -1 . (98) 


We now show that for any fixed A, N t and K. the inequality 
in ( [97| of the random quantities (mf (2)) also hold for their 
respective deterministic approximations ( fhi(z )). This can be 
argued as follows. Consider the difference 

mf(z-i) - mf (Z 2 ) = mf(zi) - fhi(z{) - mi(z 2 ) 

+ fhi(z 2 ) +m i (z 1 ) - fhi(z 2 ). (99) 

First note that since the matrix Hf has bounded spectral norm 
almost surely, it follows that, almost surely, 

lim rrii(z) > 0. (100) 

L—f 00 

From ( |100[ >, further applying the result of ( |96| in the right 
hand side of ([99]) and the inequality of ( |97| ). it follows that 

TTliOi) - ffli(z 2 ) > 0. (101) 


Thus, fhi(z) is a decreasing function of 2. We remark that the 
strict positivity of the rnf(z) in ( |100| > is essential for the strict 
positivity argument of ( |101| ). 

Monotonicity: Consider //' > fi. In this case, we denote 
m[ 1 \ 1) and fh^\ 1) as the solutions to the fixed point 
equations in ( fl2| evaluated at //' and fi respectively. As before, 
let us consider the difference between the quantities. 


7 i,j 


7 i,j 


( 102 ) 


fijiv') - _ (1) _ Wm 

We now have to show that fh[ 1 \ 1) < fnf\ 1) in order 
to prove the monotonicity result. This can be shown by 
constructing the extended matrices as in ( |93| > and noting the 
monotonicity property for the associated random quantities 
and extending the result to their deterministic approximations 
for any system dimensions. The proof is similar to the scala¬ 
bility result and hence omitted here. 


Appendix C: Feasibility Conditions 


A. Part I: Proof of Lemma 2 

Let us first consider the expression for the achieved uplink 
SINR by the ROBF algorithm given by 


A *jO) = 


|v ff h- ■ ■ I 


V fejA lyg Vi. 


>. . 112 * 


(103) 


We examine the expression for the uplink SINR in the large 
system regime. It can be shown that the expression for the 
uplink SINR converges asymptotically to the following (the 
proof of this is very similar to the convergence result of 
the downlink interference and downlink SINR provided in 
Lemma 4 and Appendix E, and hence omitted in order to 
avoid repetition): 


A i ,j (fO 


/L,J (7 




E(ra,fc)#(bj) N t P'n,kGi,n,kGi,i jTTl^ + OtiGi^jin^ 

__ (J04) 

Nt P J n,kGi,n,,kW l i + 


Let us we consider a slightly modified version of the expres¬ 
sion < n~04i ) as follows: 




E n,k jrPntkGi.ntkm'i + a,m' 


(105) 


We now examine the expression ( |105| > in detail. Recall the 
fixed point equation for the computation of m, in ( | 1 2\ . Upon 
rearranging the terms, we have 




k^n^k^i 


Nt 1 + Oi,n,kH j n,k , Wli 
n,k 


(106) 


Substituting the expressions for rrii and m! i in ( | 105| ), 

we get 

(J . . .... . (\ _lY (o-j ,rt, k Vrt,kfhi ) 2 \ 

, _. u h*,3P*,j l L Nt 2^n,k (1+cFi „ k fhi) 2 1 

m= — 7-Vy ■ ( l07 > 

' N t 2-tn,k (1 +<7i,n,kV-n.,k™-i) 2 

Multiplying and diving by to* in ( | 1Q7] >, we obtain 

a u rfi (1 -J-V \ 

V*,i,]Pi,3"h ^ N t 2^n,k (l+cr t^^ii^kfhi) 2 J 


< flQ5l > = 


N t Era, 


k (1 


Cfititjpi,j m i ^1 N t E 


N t Z-*in,k (1+cri.n.feAtn.fcTOi) 2 


1 —Ly 

1 Nt 




._r- — V *. „ 

N t Z-tn,k l+o-i, n ,kVn,kmi N t ^n,k (l+oy,„, fc /i„, fc m;) 2 

rj. • IJ- ffl (l _ — V (<7i,n,kVn,krhi) 2 \ 

\ L Nt 2^n,k (1 +<Ti, n ,kVn,kfhi) 2 ) 


,,kk'n,k'm i 




,, k /In, k 'TR'i 


1_lr _ = _ 


(108) 


We notice that that the expression at ( |108| > is exactly in the 
same form as the right hand side of the fixed point equation 
m Therefore, we can rewrite © as follows: 

— 2 

. _ __ yy- • 

ii,j 1 /A Ay — / - /nr —/ 

2^ n ,k n[ f Jj n,kGi^ n ^kGi^jm i + Gi^^m i 

Rearranging, we have 

^ ™2 _ \ ^ f^Ti,k Ay Ay _ / 

— / J ~~jTT ^ sr i,n,k^ T i,i,j'^ j i 
r )i,j , 


Gi,i,jTfli 


Vi, j 


E 7*,j f^n,k Gi^ji^G 
' 


n,k 


N t 


&i,i,jGi i.jTflf 


7 i,i,j r ^ l i 


Vi, j. 


1 i,x,j 

The equation in matrix form can be written as 

n = T(A') t n + k 

T 

where Ki = 


(109) 


and k = 


[ki,...,kjv] . The matrix T is the same as the matrix 

r defined in the description of the ROBF algorithm. The 
matrix A' is defined as follows. 


(A ') 1 * 1 ... (A') 


'M ,N 


A' = 


( 110 ) 


(A') 




(A') n,n 
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where each submatrix E C KxK is given by 


(A')S = < 


Lr- 




nn. 


j^C 1 . n 
N t ' Jr n,i,j Kjr ' 


n,n,k'm n i 


n = i, k = j 
n = i, k^j (111) 
n i. 


Notice that if the SINR targets are feasible for the uplink 
solution, the linear equations in ( |109| ) must have a solution. 
Problems with similar structure has been studied before in 
the context of power allocation in wireless networks ll32ll . It 
has been established that a feasible solution exists for such 
problems if and only if /?(T(A , ) t ) < 1 (where p is the 
spectral radius of the matrix) and the matrix I — T(A , ) t is 
invertible |32l . Now recall the matrix A defined in equation 
Observe that the matrix A and A' only differ in 
the diagonal element. Also, we can note that A < A' 
(denotes element wise inequality, refer to the notations), and 
hence TA < TA 7 . Hence, it follows from Lemma [7] that 

p(r(A) r ) < p(r(A') T ) < i. 

Now recall the linear equations for computing 8 given by 


8 = TA8 + p. 


It has been established in works before ll33l , the eigen values 
of the matrices TA and r A 7 are the same. For completeness, 
this can be argued as follows, 

|rA - ai| = |r ||a - Ar _1 | = |r|| (a - ai^ 1 ) 2 ’ | 

= |r||A T - Ar” 1 ! = |rA T - ai| = o. 


From the above discussion, we conclude that the spectral ra¬ 
dius of TA and TA 7 are the same and, hence, if p(TA) < 1 
then p(rA 7 ) < 1 and vice versa. Consequently, the matrix 
I — I A is invertible. This proves the result of Lemma [2] 


B. Part II: Feasible SINR targets for the modified system 
Consider the equation ( |25j ). Rearranging ( [25] ), we obtain 

i K , .mod™ mod 

1 _ J_ a i,i,kFi,k m i 

AT ' ^ i I _ , .mod™ mod 

iy t 7 ^ -1- + k m i 


k =1 


K 


1 


N K 


VEE 


(Tn 




mod 

k 


Nt^^l + VrnaAn)^™™' 1 
n^i 

li,k 


= — V 

N . X, 


fh? oA 


Nt k=1 1 + 7i > fc 


wXT. 


N K 


c{n)iCt™r A 


N t^ 1 £X 1 l+VrnaAn)lt™£rh™ A 
n^i 


m? oA . 


( 112 ) 


Equation ( | 1 12| > is true for all m ™ od ,i = 1 ,,N. Therefore, 
by symmetry of the fixed point equation fh™° A = rh mod 1 i = 
1,..., TV. Thus, 


(n)/OrC od = a max (n)p™£m mod = 


= («). 


'n,n,k 


~7n,k- 

(113) 


Substituting ( |113| ) in ( | 11 2[ ) yields 



+ TO,. 

(114) 


Now recall (|24|). Rearranging, we obtain 


u mod = 


7 i,j 


Vi,i,jfh™ A ' 


(115) 


The feasibility condition [Cl] requires that pf° d < oo, Vi, j 
This implies from ( |115| l that m.™ od > 0. Using this in ( |1 14] i, 
we obtain the following condition on the target SINR for 
feasibility: 


K N K 

1 y- 7 j,k \ ' 


c( n ). 


t n=l k= 1 1 + 
n^i 


7 n,k 

% ' n ' k < 1. (116) 




'Iri^k 


1) Proof of Corollary 1: Let us rewrite (| 11 4[> as 


i K 

1 1 7 i,k 

m * = 1 - AT Z. 


4 ee 


N K <Tmax(ra) 


t 1 + 7 i,k 

k—1 


T n,k 


N t ^ ^ 14- 

1 n=1 k =1 1 ' 

n^i 


»(») 


^n^k 

(117) 


The feasibility of the uplink problem in the asymptotic do¬ 
main implies that lim supjy t k^oo ^Tf < oo. This i* 1 turn 
implies that lim inf jv t ,k^oc fhf od > 0 (strictly positive). 
Thus, we consider ( |1 17| ) in the asymptotic regime, and take 
liminf^ jf^oo . This yields, 


lim inf to,- = lim inf 


i K 

, _L 7t,fc 

/V. X 


N t ,K^oo N t ,K->oo l N t X 1 + 7 ik 

k—1 

N K <Tma.*{n) 

_J_ °n,n,k > n i k 


4 Tl= 1 k= 1 1 + 

n^i 

= 1 - lim sup 'S'' —’ fc — 

N t ,K->oo LTV t ^ 1 + 7 i t k 


CTmax(n) 




J_ Vn,n,k 

AL ^ A^ 1 i Vma. 


N K 


“T n,k 


(n) 


N t ^ ^ 1 + 

n=l k= 1 1 + a nnk 
n^i 


^n^k 


(118) 


(119) 


Since lim inf .k —^oo mf oi > 0, and ( |119| l, it follows that the 
feasibility conditions are given by. 


lim sup 

>oo 


1 K 
wX 

1 fc=1 


7 i,k 

1 + 7j.fc 


1 

jv* 


JV if £V, 

EErf 

n= 1 fc = l 1 ' 

n^i 


r(n) 


T n,k 


„(n) 


7n,fe 


< 1 V*. 

( 120 ) 


C. Port ///: Proof of Lemma 3 

Let us denote /x” od = [/i™ d ,.... p,™fi A ] T and p moA = 
[/x™ od ,..., /x™ od ] T . The main idea behind this proof is the 
following: Consider the pair of vector { 7 ,/r mod } that satisfies 
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©■ We show that the when the power allocation /r mod is used 
in the original system, then 


a /i mod ?T! modl > n- ■ 


Vi,j 


where m™ dl satisfies 


1 


mf odl 


mod 


Nt 1 + <7 M ,fcM™ d ™” odl 


AT K 

EE 


_ ..mod 

^i,n,k^ nk 


l 

7U Z / Z / | I _ . .mod™i 

n=1 fe=1 1 + ^i,n,k^ n ^k m : 
n^i 


mod™ modi 


+ 


( 121 ) 


( 122 ) 


where \\ 22 ) represents the fixed point equation calculated by 
utilizing the power allocation /x mod in the original system. 
If ( | 121 1 > is true, then this implies that 7 is a feasible SINR 
target vector for the original system (since there exists a power 
allocation // rnod that can achieve this SINR target). 

We prove this as follows: Firstly, it is easy to see that 


Vmaxin)^ > CTi,n,kHZ 


since cr max (n ) > cx irl fc. Recall that the function 
fij = 7 ‘ ,J is a standard function. From ( | 123| > and the 

monotonicity property of standard functions, we have 


7 i,j 


> 


7 i,j 


o-i,ijm? od a. 




mf odl ' 


(124) 


Rearranging ( | 124[ > we have, c r J ;, ij m™ odl > (jj l 3 rri™ od ■ Multi¬ 
plying by ^™° d , it follows that 


.mod- modi 


’-.I m i 


.mod- mod 


(125) 


From ( |24| i, = 7 ij. Therefore, we can conclude 

that 


mod- modi 




> 7 i,j- 


Appendix D: Convergence of the Uplink SINR 

In this Appendix, we prove the result of Lemma[3] statement 
( |36| >. Recall that the expression of uplink SINR for the ROBF 
algorithm is given by ( |32| . We make a variable change and 
denote = *7^. This implies that 


4“ I-Nt) 


Vi,j. (126) 


Note that /i,;j is independent of the elements of the chan¬ 
nel matrix. Applying Lemma [9] to the quadratic term of 
+Iv t )" 1 x Mj - yields to. 


E 


H 


1 


+ii Vt )- 1 x Mj --— tr(s; +i w j 

. Ci w . . 

< — Vi,j 

N7 


(127) 


for k > 2, and constant Ci independent of iV t and K. 
Additionally, from the result of Lemma 6.1, l34l . we have 


E 


^tr(Sf+Ijv t ) 1 ~m l 


< %, Vi, j (128) 

N t 2 


for k > 2 constant C 2 independent of N t and K. There¬ 
fore, from fll27| ), ( |128| l and Holder’s inequality (\x + y\ k < 
2 fc_ 1 (|a:| fe + \y\ k )) we conclude that for some constant C 3 , 


E 


H 


’i 


In, ) 711., 


< 


C, 


— k J 


Vi,]. (129) 


N3 


From d 126K we have, x[ r t J CS 




T.. . . — KAvi.i) 

1 N t ) — fj. 


Moreover, at the convergence of the fixed point equation (ft 


Ti,j 


jMi, 


-. Substituting in \\29\ we have, 


E 


7i,j 


C 4 W ■ 

< Vi,j 

N 2 


(130) 


where C 4 = C^Ci i jfii j. Note that 07 ,; j is bounded. More¬ 
over for any target SINR 7 satisfying the feasibility con¬ 
ditions of limsup Nt K -toofH,j is bounded. Therefore 
lim sup iY; K ^ f00 C 4 is bounded as well. In order to prove con¬ 
vergence results, we examine the properties of the supremum 
over all the indices i,j. 





k~ 

(“) x—. 


k~ 

E 

sup 

A hjifcj) li,j 


< E E 

A hjiU’ij) li,j 



. i,3 





- 


(b) 

< NK 




(131) 


where C 5 = C\ N 8 (where 8 = jV , ' d finite value). Inequality 

(а) follows from the linearity of expectation operation and 

( б ) follows from the bound in ( | 1 30| >. Finally, we make use of 
the following inequality (which can be shown easily, details 
omitted here) 


supE [|A ij(iMj) - H,j\ k ] < E 




sup |A 




From (JT3TJ and ( |132| ) we deduce, 

C 

supE [|A - 7i,il fc ] < —fc 
ij N t 2 



(132) 


(133) 


By taking k to be sufficiently high (k > 6 in this case), the 
right hand side of ( |133| l is summable. By Markov Inequality 
((5.31) of |[35l ) and the Borel Cantelli lemma (Theorem 4.3 
of li35l ), it follows that A ?iJ (/pj) — 7 ij ———> 0, Vi,j. 


Appendix E: Convergence proof of the Downlink 
Interference 

We only focus on the convergence of inter-cell interference 
term (the convergence of the intra-cell interference follows 
in a similar manner). Throughout this section, we denote 
= Ebli EfcLi hn, 6 ,fch^ 6)fc + I. Consider the inter¬ 
cell interference term Ea l v ^fclw,j | 2 ( n 7 ^ *)■ Substituting 

for V„ jfe = n,n,k, we obtain, 

E l V ^fc h «.i.j| 2 

k 


v n,k h H ^>“ 

k 1 






( 134 ) 
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Performing a change of variable x„ ( jj = ^ph n lJ and 
x„,n,fc = we obtain, 

\ V n,k h n,i,j\ 2 

k 

= x n,i,jK ' ^n lx n,i,i- (135) 

k 

We use Lemma [ 8 ] to remove the column from the matrix 
3> n . Denoting $„ = $„- //,. / x/,', ; , we obtain 

k 

' 2 ,i,j^n {^2k ^n,k' x -n,n,k' x -n,n,kj x n,i,j 


+ x n,bf) 2 


.Ax„,., - ^tr(A) 




N t 


< 


A 


Nt' 


If 11A11 < oo, then we can conclude that 

“tr(A). 


Ax 

Vij 


N t 


N t 
(b) J_ 
N t 


H —2 

C n,n,fc^n x n,n,/c — 


X n,n,fc $ " X «,n,fe 


Let us focus on the numerator term 
% Efc tin.k (xn,n,A~ 2x n,n,k) ■ Using the fact that 


N t 


N is bounded, and applying the result of Theorem 0 eq. 
([84]), it can be shown that (the exact details of the derivation 
are omitted) 

Efc <W (*n,n,kK _ Jfk K k q n ^ k m' n 

N t ~ N t 

From the results of (| 1 36|> and (| 143|>, it can be concluded that 


■ (143) 


E l ff . 12 _ bn,kG ni ijG n} n^m' n , 

l v n,fc n ™,zj| at > n r 


N t 


(144) 


(136) 


First, it can be easily verified that the denominator term in 
( ]136| ) converges to 

(1 4” x n, i (1 4“ ^i 1 jO’ n ^i^jTfl n ) • (137) 


D. Boundedness of ||A|| 

We will complete the proof by showing that 11A11 < oo. 
First, it can be easily noticed that the matrix A is Hermitian. 
We now show that the matrix A is also positive semi-definite. 


For any vector g G C JV * xl , we examine g H Ag. 


We now focus on the terms of the numerator of ( | 136|> . Let us 
denote the matrix A = $'„ _1 <Wxn, n ,fcX^ njfc J 4>' re _1 . 
Applying Lemma |9] on the term X.U ^Ax,, ,^, we have 

Cr, 


1 ^ 5 ri,k*n,n,kX-n,n,k ^ K 's 
= ^T-g H$ n^ lx n,n,fc X ^n,feC _1 g (145) 


(138) 


(139) 


Denoting y = $' n L g, 


^2s n ,kg H yyg H = '^2s n ,k\g H y\ 2 >0 (146) 


k k 

r |2 


The proof of the fact that ||A|| is bounded is deferred till 
the end of this section (see Subsection D at the end of this 
appendix). 

Let us now focus on the term -^tr(A). First, we perform 
some straightforward manipulations. 

^(A.) = ~rtr ^ 5 n , fc x n , nifc x^ riifc ^ 

= TT^ ( (^2tn,kXn,n,kX* n ,k) <” 2 ^) 


since d rlk > 0 and g"y | 2 > 0. Therefore, the matrix A is 
semi-definite. Next, we note that (since A is Hermitian and 
positive semi-definite) 

11 A| | = A max (A). (147) 

Consider z, the eigen vector corresponding to the maximum 
eigen value of the matrix A. 

11A11 < z H $; _1 ^ 5„,fcX„,„ ifc x^ n >fc ^j $; _1 z 

= ^ ' ^rt,fc z 4 , n x„ i „ i tx n ni |.$„ z 


J^ S n,k (xn,n,kK ^n.n.fe) , (140) 


= £ 


^ n > k ) .. tf*/ t 




Mn.fcZ yi n ,u,k^n,n,k^r 


Note that j3 = max n> t jf-^ is a bounded value. Substituting, 


were (a) follows from tr(AB) = tr(BA) and ( b ) follows 
by noting tr^Aj) = X^ tr ( A i) and tr(xx ff A) = x ff Ax. 
Applying Lemma [ 8 ] to extract the column vector : ’X- n , n ,k from 
the matrix and denoting <f>" = <f>' rl - ^n,/cX„, n , fc x^ n k , 

we obtain. 


we have 


|A|| < P^2fJ, n ,kZ H K lx n,n,k*n,n,kK ^ 


r- (141) 


(1 + /^n,fc x n5n? fc^n x n,n,fc) 2 

Note that the denominator term of ( | 141 1 > converges to the 
following: 

(1 4“ dn.k^-n n k^n x n,n, k ) ^ (1 T &n,n,kB , n,k'm r i) • 

(142) 


< P l Kn,n,kX-n,n,A 1 Z 4 - /?Z H2 Z . 

k 

where the last inequality follows by noting that j3z H <&' n ~ 2 z > 
0. Therefore, 

II A|| < pZ H <S>' n ~ l (^2 Bn,k^ n ,n,k^n, n ,k + ^ 

< /3z H $' n ~ 1 z < /3A max (4>'„ _1 ) < oo. 


(148) 
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Appendix F: Convergence of the Lagrangian 

First we start with the proof of ( |40| ). Consider the La¬ 
grangian given in the form of 0 - For simplicity of notations, 
let us define 

\v H h - I 2 

IT. . A 
u l,3 — 


1ij 


4,? = E 

fc /2 


H i i 2 f(l) ^ &i,kGi,i,j(-*i,i,k'n7 j i 
I > 1 iJ — 

k^i 


r( 2 ) A V lv ff h - -I 2 f (2) - 

2-^ \ V n,k n n,i,j I > 1 i t j — 


N t 

&n,kG n jj G nn ] i Ul n 


n,k 


n,k 


N t 


n ^ i. 


E 

id 


Ai, j 
N t 


U, 


id A 1 ) r( 2 ) 


- DV - />7 - iVn 


lid 

E AiJ 
N . 


1,3 


*> 3 1,3 


U i,3 f(l) t( 2 ) 


7i,j 


7A- L ,' T\*J AT 

l-.l '■:) 


Therefore, 


lim L dl (v,/z)x lim V vf,v; . 

N t ,K^oo V ' N t ,K^oo ^ l ' J ,J 

i,3 


+E 


A 4 *,2 


AV t |_ 7 c 


^ij 7-(l) r( 2 ) at 

^ U ,2 ~ ^ iV t 


Following the argument of ( p~4] >, we note that 

TL ■' - - z£ } - iV 0 = 0 Vi, j. 


7*,j 


i,3 i,3 


as 


E v 5v*»t = EK v >„ 


^ lyg.h- I 2 

3 ju l v *,2 n *>T2l 


1,3 


1,3 


jijNt 


E 




An,k | „g u 2 

jy l V i,i* 1 *,ra,fc| 


can be shown (the details are omitted here as they follow from 
steps similar to the proof of previous results). 




°i,i,3 m i 


(1 + //,. ; rr,.,.,m,) 2 


v " 1 , . . 12 
JV, 


Aid 




1 H - 


E An,k | u 12 


(149) 


( n,k)A(id) 

E 

(n,k)^(i,j) 


l^n,k^i,i,j Oi,n,kTTli 


N t ( 1 + md a i,ijmi) 2 ( 1 + An,kVi,n,kmi)‘- 


(156) 

(157) 

:• (158) 


We examine the asymptotic convergence of the terms of the 
Lagrangian, i.e., limjv t ,A'->oo L(v, /t). By using arguments 
similar to the derivation of (|143|), the following can be shown: 


From ( |157| ) and using the fact that Jijj = a i,idAid m i’ 
follows that 


v/'.h -> 2 


^ i,z,3 11 “i 


lid N t i,J l ' l '° (1 + /ii J cr i ,j,jTO I ) 2 ” 


(159) 


(150) 


Additionally, similar to the proof of ( |150| >, it can be shown 
that 


y d N 

Z^ AT . 


1 ~ r HAid \-? r H U I 
AT, 7,2^' i,ijl 


(151) MAAd) 


An,k | - H i. 12 

jy 2 


^ ^Iv"h 


Ej 

Z^ AT 


1,3 


N t L(1 + HijCTi.ijmi) 2 (1 + fj.ijai.ijmi) 2 


(152) 


E 


fJ"n,k&i,i,j Gi,n,k’kli 


From ( |150| > and ( |152| i, we conclude that 

.I™ L DL (v,fj) x lkn (153) 

Nt,K—>oo Nt , —>• oo * ^ J 

i,j 

This completes the proof of (|40|. 

Next, we proceed to the proof of ( |4 1 [ i. Consider the La¬ 
grangian in the form of Q. Using (|5j, we can rewrite the 

term E ij V /, 2 B ' i V '- 


[n,k)jt(id] 


Nt( 1 + fiij<Ti,id m i) 2 { 1 + An,kCT i,n,krrii) 2 1 


(160) 


Therefore, we conclude that 


lim L ul (v,u)x lim V 

AC —V ATj. AA —i • ^ 


Aid N 0 


N t ,K—t 


v E 

z^ AT 


N t ,K-t oa^ N t 

1,3 

a i,i,3 m i 


'1,1, J ,,l 1 


N t L(1 + fiijaijjmi) 2 (1 + fi,.,a, , jin,) 2 


E 


(154) 


(n,k)^(i,f 


^n,k^i,i,j^i,n,k^i 

N t { 1 + flijai,ijffli) 2 (l + fJn,k(Ti,n,kffli) 2 \ 

(161) 


Substituting for = W jf-^ij in ( |154| ), we obtain Finally, we show the following result: 


E v 5 B M v ‘.2- = Efe[^.i 


1,3 i,3 

EL|v 7/ h • ■ 

At 1 


7ijivt 


{n,k)^(i,j) 


jy t 


+ E 


(1 + fJijaijjmi) 2 (1 + fiijajijmi) 2 

\ ^ l^ j n,k^i,i,j^i,n,k^'i 


(155) 


, . .. AV t (l + fiija.i,ijmi) 2 (l + fJn,kcri,n,kmi) 2 

[n,K)^{i,3) 


= o. 


We now examine the asymptotic convergence of the terms on 
the right hand side of (| 1 55|>. The following convergence results 


(162) 


The proof of this result is provided next. Consider the sum of 
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the first and last terms of (|162|i. 


E 




n,k 


N t { 1 + / 2 jjCr iii jTO i ) 2 ( 1 + / 2 „,fcCr.i,„,feTOi ) 2 


+ 


" v i 


' i,i,j ,,b % 


(1 4- Hij(Ti,ijmi) 2 


1 ^ r f-^n,k^i,n,k 




(1 + mj(7i t ijmi) 2 \ N t f^{l + Hn,kCTi,n,kmi) 

(163) 

Multiplying and dividing the right hand side of ( | 163| > by m 7 , 
we obtain 


' i,i,j nh i 


nrii{ 1 + 


1 f-ln,k&i,n,k‘Wli 

N t nk ( 1 + 


(164) 

Recall the fixed point equation for the computation of rhi in 
©• Upon rearranging the terms 


m, 




Oi ,n ,k Mn , k^i 




(165) 


Using rrii for the terms inside the brackets of the right hand 
side of (| 164|>, we have. 


— Y- 

N t *-! 1 


& i,n , kl^n^k^i 


— y 

N , 


T ( ^i,n 1 k^n,k’kl r li Nf. (1 T H"n,,kGi,n,k'kkli) 


m-i 2 


1 Ai,n,k^n,k m i) 

N t nk (1 + (J i,n,kHn,kfh i )‘ 1 ' 


Rearranging the expression for m' in 

= ,_L 

fh ’. N , 


& i,n, k^-n.k^i 

N t V 1 + Vi,n,kHn,krhi 


we have, 
2 


Using ( | 166| > and ( |167| i in ( | 1 64[ >, it can be verified that 
( fl64| = 


(1 + cn , n , kHn , km i )' 1 ' 

From ( | 168| >, the result of ( | 162| > follows. 

Finally, using ( |162| > in ( |161| ) we conclude that 

IH,jN 0 

~nT' 


lim L VL (v,n) 

N t ,K^oo N t ,K-to 


lim > 

t ,K—i^ ‘ ^ 

This completes the proof of (|4T|>. 
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Appendix G: Convergence of uplink power 

ALLOCATION 

The proof proceeds by an approach similar to fm of 
finding the upper and lower bounds of the optimal solution 
and show that these bounds coincides asymptotically with 
the solution obtained by the ROBF algorithm. Recall the 
fixed point equation for the computation of the uplink power 
allocation in the CBF algorithm in (pj. It can be rewritten as 

fiM = t ■ g -Vi, j (170) 


7»„ 


'sX.aY • Iv,) 1 1 •.../ 


Also recall from the arguments of a that /,; 7 (/3) is a standard 
function. 

Let us denote fki,j{Sij), as the solution provided by the 
ROBF algorithm with target SINR 7 ij + Sij Vi, j. where 
Sj j > 0 is a small positive constant (i.e. the solution pro¬ 
vided by the fixed point equation ( [TT| with J replaced by 
7 ij + Sij , Vi, j). We now examine the achieved SINR in the 
uplink with a power allocation of /ii,j(8i,j), i.e.. 


A. 


1,3 (m(^)) 


N t 


h H - ■ 


(£f 




tljvJ-'h. 




(171) 


Similar to the result of Theorem [3] (convergence of the uplink 
SINR), it can be proved that 


A ,,(jaj (AjY) - ( 7 i,j + Sij) - Y - > 0 , Vi, j. (172) 

Nt,K^oo 

Since Sij > 0, we have 


■Aij (aLj — d? Vi,j 


(173) 


Since the achieved SINR with an uplink power fa j(6i j) allo¬ 
cation is asymptotically greater than 7 ,j, the power allocation 
IM,j(Si.j) is a feasible solution to the uplink problem with 
target SINR 7 ij. From the monotonically property of the 
standard function /(/3), we can conclude that 


> \,j Vi,j. (174) 


Similarly, let us define Sij) as the solution to the 

fixed point equation ( fTT] > with target SINR 7 — 8ij, ^i,j- 
By similar argument as above, we have 


A ( li,j 




a . s . 

N t ,K^oo 


■» 0, Vi,j. 


(175) 


Using the definition of almost sure convergence, we can write 
that for Cij > 0, there exists N' such that for all N t > N', 
we can have 


A*,'/ {jJ-i.j ( I ^ (176) 


If we can select a value of G,i < Sij, then there exists a large 
enough N" such that for all N t > N", the following 


Nt 


i H - ■ 
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(Sf 




+ IjvJ 
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h, _ j 7jj, \/x, j 

(177) 


holds true. Equation ( |177| l implies that 
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hf 2 . 7 (Sf 


= fiAvi- 6 ))- 


s(-« 


nj N t 


iNj-'h 
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v«, j. 

(178) 


In other words, an infeasible point for the standard 

function 

Let us consider that we want to find a fixed point to 


Pi,j = Vi, j. (179) 

The fixed point corresponding to this is the optimal uplink 
power allocation of the CBF algorithm A,;. r We start with 
an initial value of A 3 = /u ; 7 (—Aj)- Let us consider the 
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iterations of = f t j((3*'). with /3° = Rewriting 

the hrst iteration, we have. 


(a) 


V«,i, 


\ “■ 

Mi, 7 / 'i ,7 

J 5,7 N t ,K^-oo 

Similar to the proofs as before, 

Ei ,j SiJ 'Vj o.s 


> 0 Vz, j. 


N t 


N t 


N t ,K-¥ oo 


>0 Vi, J. 


Appendix H: Impact of Imperfect CSI and Pilot 
Contamination 


Recall the uplink receive filer with CSI estimate in (65 
First note that for the MMSE estimate of the form (62), 

hj,n,fcj n/i and hi.i,k are related as 


^■i,n,k 


&i,n,k j 


H,i,k 


Using (|185) in ((66), we have 


N K 


— ^ ' 'y ] ^ + I 


n— 1 k—1 


= £ 




= i,fc x M,fc + 1 
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k= 1 
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where 
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l tVf 


V est h ■ — 
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— ,/A v[/ ^ x - 

— \f Ul >J ^ i 


First note that using Lemma [8] we can remove the column 


Xjj.j from the matrix d/ 


-l 


(180) 


as follows: 


where (a) follows from ( | 178) . Also note that this is true for all 
i,j. Using the fact that the function / is standard and hence 
monotonic, the sequence j3j j, t = 0. 1,2 ,..., monotonically 
increases and converges to X l:J . Therefore, it follows that 

^i,j • (181) 

Therefore, from the above arguments, we conclude that 

a.s. a.s. 

X/.,j ) 5) l-Li-j ) Vi,j. (182) 

By taking A,;arbitrarily small, we have 


( p~89) = 


1 5ij'kf i J T' i 


1 + C'„/ 


j 


i+Cri,,j»'r 


(190) 


where m® st can be computed as in (j69). The result in ( |190) 
can be derived as follows: First, we look at the numerator term 


X V. . U/ ^ v - 
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(183) 
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(191) 


where (a) follows from the results of Lemma [9] Theorem [5] 
and noting that the column vectors x, ^- V & / t, n-r r are 
independent of Step (6) follows using the relation in 

((64). Similarly, the denominator term converges to 


1 + €i,jX^i j'&i 1 ’ 


j 


1 + 


We now investigate the asymptotic equivalent of the sum of 
interference power at UTjj from the BS of cell n, given by 


(185) 


J2 l v n‘ft ffh n,i,il 


°"> fc uH 


7 j jy 2 n ,n,k^ n 


\H 


US' 1 h„ 


— y J ^n,k x n,i,j^n X n>nj fcX nJl x n,i,j- 


We note that the CSI estimate of UT ij is contaminated from 
the pilot signal of UT n j, (n / i). Therefore, we first analyze 
the interference signal from UT n j, (n / i). as follows: 


(186) 


(187) 
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Additionally, let us denote 

K 

+ 1 N I . 

k=l ,k^j 

K 

k=l,k^j,k^l 

First, we focus on the asymptotic equivalent of the useful 
signal term. 


(1 + £m,, 


(192) 
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where in (a), we have removed the column x w m 
matrix 'L“ 1 using Lemma [8j and the result of ( b ) is derived 
similar to ( | 190) . 

Next, we consider the sum of interference signals from rest 
of the UTs in cell n, i.e.. 


(188) 


E S n,kXn,i,j n 1 X n , iii 

= X„ 
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where B nj - = Yjk^j S n,k^-n,n,k^n, n ,k- Usin s Lemma 
can decompose ’L“ 1 as 
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ra n ’J nt - i-H 
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Using ( |194| > in ( |193| ), we obtain 
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(197) 


Consider the term ( ]195| ). It can be proved that the spectral 
norm of is bounded (following the steps in 

( |148| ». Therefore, using 

'tr(^.B„,^-i). (198) 


X ra,*J^nj P raj^ra,j X ra,ij 


N t 


Using some straightforward steps (similar to the derivation of 
(|140|)), it can be shown that 
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Nt 




i 

Nt A"' 


'If., 2 x„ 


n,n,K n,j 




The right hand side of ( | 199[ > can be further analyzed 
follows: 


(199) 
as 
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( 200 ) 


where ( m/ n ) est can be evaluated as in and the result 

of ( |200| ) has been derived following the approach of ( | 1 43[ >. 
Similarly, it can be shown that 
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Combining the convergence results of the useful signal and the 
interference signal terms derived in this Appendix, we obtain 
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